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ABSTRACT 


One method of river closure is by end-dumping loose non- 
cohesive materials from one or both banks progressively into the flow. 
This study was undertaken to experimentally develop comprehensive 
hydraulic design procedures for end-dump closure of alluvial channels. 

A numerical measure defining stability, termed "efficiency 
of closure material", was introduced to denote the ratio of dumped 
material remaining within the dump zone at any stage of closure to 
the amount of material dumped. Efficiency of closure material, scour 
depth and backwater rise were correlated with the variables describing 
the characteristics of the fluid, the flow, the bed, and the closure 
material. 

The most significant non-dimensional variables were found to 
be the Froude number of approach flow and the contraction ratio. 

Experimental data was obtained from over 120 test runs with 
Froude numbers of 0.50, 0.29 and 0.10,using sand of median size 1.2 mm, 
0.60 mm and 0.25 mm as bed material. Pebbles and rock chippings of 
6.60 mm, 18 mm and 25 mm were used as closure material in a 7.5 foot 
wide flume fifty feet long. 

Neither mean velocity through the constriction, nor velocity 
against the boundary was found to be a practical criteria for stability 


of the closure material. The stable size was found to correlate with 
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boundary shear stress as determined from measured velocity profiles. 
Design curves for stable material size, maximum scour depth, 
backwater rise and rate of scour are presented. The design curve 
presented for stable material size is based upon a closure material 
efficiency of 90%. Comparison with the few available prototype 
closure data and with other studies of scour at constrictions are 


given. 
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Ce = generalised Chezy coefficient for seepage (Eq. B.6) 


C. = discharge coefficient for constriction on mobile bed 

C, = Chezy coefficient in open-channel flow 

CC, = coefficients in "velocity-stable size" relationships 
(Eqs. 2.4 and 2.8 ) due to Izbash 

D = representative size of closure material = D556 for sizes 


above 5 mm 
Deg 2D75 D356 D1¢ = grain sizes of which the given percent by weight 


of the closure material is finer 


D, = stable size of closure material from critical velocity 
formulae 

Dd. = stable size of closure material from shear stress con- 
sidering drag and lift 

Dy = stable size of closure material from shear stress con- 


Sidering drag only 

d = representative size of bed material = deo for sizes up to 
2 mm 

dog deg 2dg5 2d) ¢ = grain sizes of which the given percent by weight 


of the bed material is finer 


d. = limiting depth of scour below original bed level 

dy odta = depth of scour at time t during active phase of scour (Fig. 6.11) 
e = base of Naperian logarithms; porosity of coarse fill 

= = Froude number of normal (approach) flow = V/yYgh 

Fy = drag force on a stone (closure material) 
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J, dood and Jy = coefficients in relationships of closure material 


stability in eqs. 3.2 through 3.15 


J} +d, = exponents in sical - m" relationship (eq. 3.6) 

5354 = exponents in sii - m" relationship (eq. 3.26) 

K = a coefficient for closure material in Sandover's logarithmic 
velocity formula (eq. 2.16) 

Ke = coefficient of turbulent infiltration through closure 
material 

Kn = a multiplying factor in scour relationship (eq. 2.21) 
due to Izzard and Bradley, and Ahmad 

K. = Strickler's coefficient 

Ke = surface area per unit volume of solids (eq. B.2 ) due to 
Wilkins 

kK. = Nikuradse effective sand grain roughness 

E = average width of closure dam along flow 

Ly ots = width of closure dam at base and at intermediate level 

Loy = length of backwater reach 

Li5 = length from closure dam to the maximum backwater section 

L 53 = length of contraction of live stream from dump line to 
vena-contracta 

Lo, = length of expansion of live stream beyond vena-contracta 

M = opening ratio = 1 - m 

M. = actual or representative opening ratio (Biery and Delleur's 


equation 2.39) 
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contraction ratio = (B-b)/B 

actual contraction ratio 

Manning's roughness coefficient 

an exponent in Garde's scour equation 2.27 

normal (approach) discharge 

seepage discharge through closure dam = Q - On 

measured discharge through the constriction over channel 
width b 

unit normal discharge 


unit seepage discharge through closure dam = 1, eh Be: 


= unit seepage in trapezoidal part and triangular part 


of closure dam (eqs. B.9 and B.10) 
unit sediment discharge for normal flow 
unit sediment dischrage at section 33 (vena-contracta) 


out of scour hole for constricted flow 


1249243 and q4 = unit discharges at sections 11, 22, 33 and 44 respectively 


Reynolds number (Vh/v or VD/v or Vd/v) 

hydraulic radius of the voids due to Taylor (eq. B.2) 
hydraulic mean radius for normal flow 

boundary Reynolds number = u,d/v 

ratio of depth of local scour at a constriction to depth of 


scour in equivalent long contraction (Laursen) 


Py oNo sks and nae coefficients in turbulent seepage eq. B.] 
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slope of the end-dump face = coto 
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Vy 


Vy Vo .V3 and V 


coefficients in time-history of scour relationships (eq. 
6.2 ) of Govindrao and Sharma 

time of scour 

characteristic time for scour to attain practical 
equilibrium 

time from commencement of scour to attain 90% of limiting 
scoured depth (h + d.) 

free stream velocity in the flow direction 

velocity at depth y above theoretical bed 

velocities at depths Y) and Yo for boundary layer flow 
velocity measured at 0.35 kK. above theoretical bed 

shear velocity (=/ Teipe) 

mean velocity of normal (approach) flow 

mean seepage velocity through porous fill 

velocity over closure material on a transverse-dump 

dam 

= minimum and maximum values of Vy that a stone can with- 
stand (eqs. 2.4 and 2.8) 

mean seepage velocity through voids in the porous fill 
are velocities at section 1-1, 2-2, 3-3and 4-4 (subscript 
e refers to velocities after equilibrium scour is attained) 
velocity at section 3-3(vena-contracta, location of deepest 
scour) at time, t 

Shear velocity 


critical shear velocity for bed material 
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Shear velocity at time t 

Shear velocity at section 3-3 (scour hole) after equilibrium 
shear velocity at section 3-3, at time t, during active 
scour phase 

maximum velocity occurring against the end-dump face 
submerged weight of stone (closure material) 

fall velocity 


linear dimensions along scour hole 


F eCiareaccertocic lengeun.O0 scour nole at, cime ti(hig., onl) 


= depths of flow normal to the boundary at which the 
instantaneous velocities are U, Uy, and Up 

Cartesian co-ordinates for live stream contraction (Fig. 5.9) 
Cartesian co-ordinates for live stream expansion (Fig. 5.9) 
an expression relating unit seepage discharge Ve to 
coefficient of turbulent infiltration K- (eq. B.11) 

a coefficient relating qo and 5, in sediment-transporting 


flow (eq. 3.34) 


angle of repose of closure material in water 
energy correction coefficient 

momentum correction coefficient 

unit weight of closure material 

unit weight of bed material 


unit weight of water 
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efficiency of closure material 

an exponent in turbulent seepage equation B.1 

kinematic viscosity 

specific density of closure material 

specific density of bed material 

specific density of water 

submerged specific density of bed sediment = Penk ey 
geometric standard deviation of particle size (subscript 


d refers to bed material, subscript D refers to closure 


material) 


intensity of shear stress 


critical shear stress for granular material on horizontal 
bed 

critical shear stress for granular material on a sloping 
bank 

local intensity of shear stress on channel periphery or 
on end-dump face 

maximum intensity of shear stress on end-dump face 
maximum intensity of shear stress on bed of the constricted 
channel 

average intensity of shear stress on uncontracted channel 
periphery 

intensity of shear stress on bed at time t from commence- 


ment of scour 
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= Shields dimensionless shear stress = Tey Noweoy Jc or 
To /(op-0,,)D 

= angle of inclination of end-dump face with horizontal 

= angle of inclination of the downstream slope of a 
transverse-dump dam (Fig. 2.1) 

= angle of inclination of Fe with horizontal (Figuuce |) 

= angle of inclination of Fy with horizontal (Fig. 2.1) 

= angle of repose of material | 

= overall angle of contraction and expansion of the 
live-stream boundary with flow direction 

= local angle of contraction at the dump line 

= local angle of expansion just downstream of the vena-contracta 

= total volume of material dumped over area 12341 (Fig. 3.1) 

= useful volume of material remaining in closure area 
IZ4 Ei Geese) 

= empirical coefficient in backwater rise equation 2.3] 


for rigid-bed channels, due to Liu, Bradley and Plate (1957) 
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CHAPTER I 
INTRODUCTION 


1.1 Closure Methods 

Prior to the construction of a major hydraulic structure on 
a flowing river a closure dam is usually built across the river, to 
divert the flow through an alternate waterway and dewater the work 
area. Closure dams are generally built with coarse, granular material 
varying in size from pebbles to loose rockfill. However, the use of 
earth or sand is possible where the situation permits. Closure of a 
river is accomplished by one or a combination of the following two 
basic methods: 

(1) Transverse-dump method: - Closure material is dumped in 
more or less successive lifts uniformly across the river, until the 
closure dam emerges above water. (Fig. 1.1a). 

(2) End-dump method: - Closure material is pushed progressively 
from one or both banks to ultimately block the channel. (Fig. 1.1b). 

A partial closure by the end-dump method is often used for 
building cofferdams for auxiliary works such as a powerhouse, a spillway, 
or a bridge abutment close to one bank. Construction of a spur ona 
flowing river is another example of partial closure by end-dumping. 
Diversion is not provided in such cases. 


The construction engineer is required to assess in advance: 
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RIVER CROSS-SECTION 
AT CLOSURE SITE SECTION 1-1 


(a) TRANSVERSE DUMP 


Contraction Ratio m= ae 


Opening Ratio M= = 


PLAN OF RIVER AT 
CLOSURES SITE 


SECTION 2-2 


Note: I, I, IT represent stage of 
closure 


(b) END DUMP 


FIG. 1.1 DIAGRAMMATIC ILLUSTRATION OF CLOSURE METHODS 
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(1) the size of stable material at every stage of closure or 
alternatively, given the size of available closure material, to assess 
the critical height of the dam for transverse dump type or the critical 
gap width for the end-dump type, 

(2) the quantity of material for the entire closure dam. 

In the case of end-dump closure of alluvial rivers, consider- 
able bed scour occurs at the dump area and downstream, particularly at 
high constriction. In such cases, the determination of the quantity 
of closure material becomes more involved. 

Design procedures to arrive at the stable size of closure 
material at any stage of closure for both transverse and end-dump dams 
have until now (1972), been usually based upon the critical-velocity 
concept, the sequence being, 

(1) assess the maximum velocity that occurs over the body of 
a transverse-dump dam or through the opening in an end-dump dam. 

(2) compute the competent size of closure material on the basis 


of this velocity. 


1.2 Characteristics of Closure Methods 
1.2.1 Transverse-dump Method 

From a trestle bridge or a cableway or a barge system, transverse- 
dumping is proceeded with uniformly across the full width of the river. 
The material upon dumping assumes and maintains a more or less trapezoidal 
shape under water (Fig. 1.la). Continued dumping results in flattening 


of the downstream slope due to increased differential head and increased 
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4 
flow velocity over the dam. The pertinent feature, however, is that the 
flow maintains a more or less uniform depth hy» and a uniform transverse 
velocity, Vi» over the closure dam. The stability of the closure 
material is then predicted by correlation with the flow velocity close 
to the particle based on the studies of Blanchet (1946,1947), Izbash 
(1936), Izbash and Khaldre (1970), Izbash and Lebedev (1961), Olivier 
(1967) and design charts of U.S. Corps of Engineers (1970), all of which 
derive from the familiar Brahms-Airy Sixth power law (Leliavsky, 1955) 
for flow on a flat bed. 

Straub (1953), and Neill (1967) have improved the procedure 
by relating the stable size on a channel bed to the depth of flow as 
well as to the mean flow velocity. It is indisputable that the mean 
velocity in a deep flow has to be greater than in a shallow flow to 
displace the same size of particle, because movement depends pri- 
marily on the shear or tractive stress on the stream bed. The mechanics 
of movement has been well summarized by the ASCE Task Committee on 
Sedimentation (1966) 

"Thus it is seen that mean velocity alone cannot express the 
scouring action of the water at the bed and that to com- 
pletely specify conditions the depth must also be given ... 
The advantage of using shear stress to specify critical 
conditions is that the one quantity suffices whereas if 


velocity is used one must also report the depth or the 
position at which the velocity is observed". 
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1.2.2 End-dump Method 
(a) Closure of rigid-bed channels 

Progressive lateral dumping to constrict the river from one 
or either banks dispenses with a bridge or a cableway and is a simple, 
economical scheme for closure of wide rivers (Fig. 1.1b). The definition 
sketches shown in Fig. 1.2 indicate the flow pattern through an end- 
dump constriction. These sketches illustrate the difference in water 
surface and bed configurations for closures on rigid-bed channels and 
alluvial channels. Increase in contraction ratio m, results in in- 
creased differential head and increased flow velocity through the con- 
striction. 

Lane (1919-20), Kindsvater and Carter (1955), Valentine (1958), 
Sandover (1969), and Izbash and Khaldre (1970), studying flow through 
open channel constrictions on rigid beds, have developed empirical 
coefficients for discharge, head loss and backwater rise. Using these 
coefficients mean flow velocity through the constriction (section 2-2 
in Fig. 1.2b) can be calculated approximately, but no information was 
available for the velocity distribution either through the opening 
or against the advancing boundary (shown shaded in Fig. 1.2a) as the 
flow negotiates the constriction, Izbash and Lebedev (1961), and 
Izbash and Khaldre (1970), however, suggest that flow concentration 
to the extent of 20% over the mean velocity occurs against the end- 
dump face. They recommend the closure material to be competent against 


the increased velocity (using empirical correlations of transverse- 
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O  sLoPE ie ] h 5 3 Normal water Surface 
Water Surface for rigid bed 


Water Surface for mobile bed 


depth-h 


Normal Flow 
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FIG. 1.2 DEFINITION DIAGRAMS FOR END-DUMP CLOSURE 
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dump dams). For the curvilinear, rapidly accelerating flow through an 
end-dump poneerieri en: this procedure is of questionable accuracy. 

Pariset and Hausser (1959), and Sandover (1970) are probably 
the only researchers to have specifically reported on the stability of 
rockfill in end-dump closures on rigid-bed channels. They have pre- 
sented empirical formulae for determining the stable size from relevant 
flow parameters and not by direct correlation with the flow velocity. 

Pariset and Hausser's procedure involves predetermination 
of the head loss and the maximum backwater rise at the constriction. 
They have not suggested any method of evaluating these and presumably 
empirical coefficients of other researchers would have to be applied. 
Sandover's procedure involves a trial and error solution. As such 
both the procedures have limitations in practical application. 

(b) Closure of alluvial channels 

The coefficients of discharge and backwater determined for 
a constriction on a rigid-bed cannot be used to determine the velocity 
through an end-dump constriction on a mobile bed because of the 
possibility of bed scour. Sandover (1969) has determined similar co- 
efficients for constrictions on a sand-bed channel in a laboratory 
flume (2 ft. wide), but his studies being inconclusive, the coefficients 


cannot be considered generally applicable. 


1.3 Purpose and Scope 


From the foregoing discussions it is obvious that knowledge 


even of end-dump closure on rigid-bed channels is by no means complete 
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and no specific procedure for alluvial channels has so far been form- 
ulated. 

In the majority of cases where river closures are required, 
the bed is alluvium underlain by bedrock. The most general case 
is that in which the bedrock lies below the deepest scour level re- 
sulting from the constriction. This study was therefore framed to 
deal with the hydraulics of end-dump closure of alluvial channels. 
The specific objectives of this study are: 

(1) to study the characteristics of flow through end-dump con- 
strictions on an alluvial bed and to develop design curves which will] 
directly determine the required stable size of closure material at 
each stage of the closure; and 

(2) to develop design curves in respect of other associated 
phenomena of interest to the design and construction engineers, such 
as the configuration of the livestream boundary, the extent, pattern 
and time-rate of bed-scour at the dump line and the maximum backwater 
rise. 

The key fact, however, is that all these flow phenomena 
mutually interact and the stability of the closure material is a com- 
posite effect. It needs to be mentioned that the mechanics of sedi- 
ment movement in uniform flow over an alluvial bed is not fully under- 
stood yet; the difficulties and limitations inherent in the analysis 
of curvilinear, accelerating flow resulting from end-dump constrictions 


on alluvial channels, are thus obvious. 
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1.4 OQutline of Contents 

Previous studies relevant to end-dump river closures are re- 
viewed in Chapter II. Studies on material stability in transverse- 
dump closures by Izbash (1936,1970) and the extension of such studies 
to end-dump closures by Sandover (1971) are examined. Specific studies 
either on bed scour or on backwater at end-dump constrictions, except 
for some qualitative informations by Sandover (1969,1970), do not 
appear to have been reported. Scour at an end-dump constriction is a 
combination of local and general scour. Investigations on local scour, 
such as due to bridge abutments or spur-dikes and on general scour due 
to channel contractions are reviewed. Studies on backwater due to 
different constriction geometry both on rigid-bed and alluvial channels 
are reviewed. The literature review revealed the lack of complete 
design procedures for end-dump closures, although valuable contribu- 
tions have been made notably by Izbash and Sandover. 

In Chapter III, comprehensive analyses of material stability, 
scour and backwater are made considering all the pertinent variables 
involved in an end-dump closure, such as the properties of the flow, 
fluid, closure material and bed material, and the constriction geometry. 
The role of shear stress in mobilizing the closure material] and in 
scouring the bed, is specifically considered. Dimensional analysis 
is applied to derive non-dimensional parameters in respect of stability, 
scour and backwater. These analyses led to the formulation of the 


experimental programme for determining the parametric relations. 
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The experimental equipment, materials and procedure are de- 
tailed in Chapter IV. The tests were conducted in a 7.5 ft. wide 
flume and comprised over 120 runs utilising three sizes of closure 
material of 0.021, 0.059 and 0.082 ft., and three sizes of sand of 
1.2 mm, 0.60 mm and 0.25 mm as bed material. The normal flow over 
the alluvial bed transported sediment in certain runs. The bed forms 
and the characteristics of normal test runs are described in Appendix A. 

In Chapter V, the characteristics of flow through an end-dump 
constriction and its inter-relation with the stability of closure ma- 
terial are examined. Based on the analysis of data, design curves 
in respect of the flow contraction and expansion, and material stability 
are presented correlating the essential non-dimensional parameters 
developed in Chapter III. The hydrodynamic effect of the flow is 
analysed, giving due consideration to the reduction of the discharge 
through the constriction on account of the turbulent seepage through 
the porous closure material. An analysis of the turbulent seepage 
in the light of a non-linear flow law is made in Appendix B. Although 
correlation of the experimentally observed stable size with the velocity 
against the closure material was found to be unrealistic, a similar 
correlation with the experimentally observed shear stress against the 
closure material as shown in Appendix C, yields realistic results. 
The applicability of the design curves on stability developed from 
studies using rectangular and trapezoidal channels in the laboratory 


is tested against the few actual river closures for which data is 
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available. 

In Chapter VI, analysis of the data in respect of scour and 
backwater, leading to the design curves is presented. The results of 
some check runs with coal and pebbles as bed material are discussed, 
indicating the satisfactory agreement of the experimental data with 
the analysis. Results of other researchers are compared with the 
findings of this study and satisfactory agreement is noted. A short 
summary, conclusion and recommendation for further studies is contained 
in Chapter VII. The complete experimental data is detailed in a tabu- 


lar form in Appendix D. 
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CHAPTER II 
REVIEW OF LITERATURE 


2.1 Stability of Closure Material 
2.1.1 General 


In order to be stable on the body of a closure dam the closure 
material has to withstand the hydrodynamic forces of the flow. Studies 
concerning the stability of rockfill forming a transverse-dump dam are 
based on the hydrodynamic effect of uniform or gradually varied flow on 
a level bed of coarse, non-cohesive, granular particles. Such analyses 
have been generally extended to end-dump closures. The available liter- 


ature pertains to closures on rigid bed channels only. 


2.1.2 Analysis for Transverse-Dump Dams 
Izbash (1936) and Blanchet (1947) are probably the earliest in- 


vestigators to have analysed the stability of stones dumped to success- 
fully close a channel transversely. Both considered the flow velocity 
over the closure material to be the most significant parameter. Their 
procedures and final equations correlating velocity with stable size 
are fairly similar. 

Following Izbash (1970), the stability of a stone experiencing 
a velocity Vy (Fig. 2.1 ) requires analysis in two possible configura- 
tions on the body of the fill. For simplicity an equivalent spherical 
diameter D, is considered as the representative linear dimension of the 


Stone. 
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STONE 1 


We 


(a) STONE ON CREST OF FILL 


(b) STONE LYING PROTECTED ON SLOPE 


FIG. 2.1 FORCES ON A STONE ON THE BODY OF A TRANSVERSE-DUMP DAM 
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14 
In the first configuration (Fig. 2.la) the stone 1, is lying 
exposed on the crest of the fill and is liable to slide down or over- 
turn about an edge by the drag of the flow. The forces acting on the 
stone are 
1) Drag, | 9 
nD° t 
aye) sea) (2a) 


in which Cy is a drag coefficient which depends upon the shape of the 
stone and a characteristic Reynold's number V,D/v; v being the kinematic 
viscosity of water; Yue the unit weight of water. 


2) Submerged weight of the stone, 


n° 
We= (ee) (yp-¥y) (2.2) 


in which Vp = unit weight of the stone. 


3) Frictional Component of Wo» 
Pree p(1>) wares (2.3) 
f 6? \Yo Vw 


in which f is the coefficient of friction of stone in water. 
From statics, both for the limiting condition of equilibrium 


either by sliding or by overturning, algebraic manipulation will yield 


Vv. = f2(2-4] (2.4) 


t represents the minimum flow velocity - capable of dis- 
min 
placing the stone and C, is a coefficient which is a function of Cy and 


in which V 
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f. Izbash evaluated C, empirically as 0.86. 
In the second configuration (Fig. 2.1b) the stone 2, is lying 


in a pocket after sliding down the slope of inclination 6, . Vi. is 


max 
referred to the maximum stream velocity that this stone lying on the 
Slope surrounded by other stones, can resist until limiting equilibrium 
is attained. The forces acting on the stone are 


1) Drag acting at an angle 6, to the horizontal, 


: 2 
2 V 
ce 7D t 
2) Submerged weight of the stone, 
We = 2 (yp-y,,) (2.6) 
Sieecca Daw 5; 
3) Frictional component of Wo > 
F = nD? -y ). Cos 6 (277) 
i Biel ‘we 2 ; 


in which Fe acts at an angle 85 to the horizontal. 
Again from statics, for the limiting condition of equilibrium 
either by sliding or by overturning, algebraic manipulation will yield 


S01 * (2.8) 


in which Cy is a coefficient obviously higher than C, and depends upon 
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Cys £5 81> 8, and 8. and has been evaluated experimentally by Izbash to 
beg1220? 

Many rivers have been successfully closed by transverse 
dumping both in the U.S.S.R. and U.S.A. applying Izbash formula 
(Torpen, 1956, Linford, 1967). 

Straub (1953) analysed the material stability from the theory 
of transportation in a wide channel. He evaluated the shear stress 
on the body of the dam for uniform flow and equated this to the critical 
tractive stress for the closure material. Referring to Fig. 2.1 


shear stress on the body of the dam is, 
CY hy mei (2.9) 


in which i = slope of the water surface. 
Tos critical shear stress (Shields, 1936) for the closure material 


is given by 
ay 0.06 (Yp-Yw? a (2.10) 


The slope of the water surface can be expressed in terms of 


flow velocity by applying Manning's formula, 


y= 1986, 2/3 | gle (2.11) 


in which n = the coefficient of roughness and is given by Chang (1939) as 


n = 0.0432 p'/é (2.12) 
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Expressing equation 2-9 explicitly in terms of Yw? ne, Ve and: Deand, 


Cheat & 
algebraically manipulating Vis the mean velocity at which the closure 
material would get entrained is given by 
We ] 
Yo Yw ; ne ie 
tbs = 8.45 Eceeaperd (—) ae) k2ahe) 


" D 
In a transverse-dump dam the flow depth hy and velocity ve 
can be reliably computed, treating the flow as a friction controlled 


varied flow, from the energy and momentum principles. 


2.1.3 Analysis for End-Dump Dams 


In a transverse-dump dam, the closure material lies on a 
level bed, over which transversely uniform flow occurs. In contrast, 
in an end-dump dam, curvilinear accelerating flow occurs against the 
material lying on a slope (Fig. 1.2a). In spite of this, Izbash and 
Lebedev (1961) recommend the coefficients, correlating the stable size 
with velocity developed for transverse-dump dams, be applied to end- 
dump dams, although no information of laboratory studies by them, justi- 
fying such a procedure, are readily available. They however, recommend 
an arbitrary increase of 20% over the mean velocity through the opening, 
which accounts for flow concentration at the end-dump face, to be ap- 
plied for predicting the stable size. 

Izbash et al (1970) have presented non-dimensional curves for 
computing the discharge coefficient Cas maximum backwater rise h*, and head 
loss at the contraction along the centre line of flow AH for open 


channel constrictions (Fig. 1.2b). The mean velocity of flow through 
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the contraction which is to be closed by rockfill is then given by 


1/2 
qo  Cglég H] 
eee (2.14) 


2 H-AH 1-AH/H 
in which H, the depth of flow just upstream of the contraction is h + h*, 
and qo is the unit discharge through the opening. 

Equation 2.14 can be obtained by assuming weir flow through 
the contraction. The stable size corresponding to 1.20V., can then be 
obtained from equation 2.4. For determining the stable size in an 
end-dump closure from Straub's equation 2.13, a similar semi-empirical 
approach as suggested by Izbash for obtaining the velocity and depth of 
flow will have to be resorted to. Depth of flow on the sloping boundary 
being variable, determination of the most critical combination of depth 
and velocity will involve a trial and error approach. 

Sandover (1971) has recently developed two specific design 
procedures to arrive at the critical gap width for a particular size of 
closure material. In the first analysis,termed critical velocity method, 
Sandover assumed the instability to be caused by mean velocity through the 
gap. However, he experimentally determined the constant in equation 2.4 
for computing the stable size instead of using Izbash's value of 0.86. 
For determining the mean velocity through the gap, he verified and 
accepted Valentine's (1958) discharge coefficient Cy for vertical sharp 
edged slot in rectangular channel as applicable to flow through an 
end-dump dam. The coefficient Cr» correlating the experimentally observed 


stable size with the mean velocity has been determined and curves pre- 
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sented by him. He has formulated a governing equation for stability to 


be balanced, in the form 


9 Cy a0) er 


6S 
2g Db“! 3 


I (2m) 
Pp-Pw) 
in which C, is Sandover's empirical coefficient, Q = the discharge, 
p = specific gravity; suffix D refers to the dumped material and 
suffix w refers to water, D = Do. = the nominal diameter of the dumped 
material defined herein as the sieve size passing 75% of the sample by 
weignt and b = the gap width. The solution does not require determin- 
ation of the flow velocity, but involves a trial and error procedure. 
His second procedure termed the logarithmic velocity method 
iS an improvement based on the concept that material entrainment starts 
after the maximum shear stress attained on the advancing boundary ex- 
ceeds a critical value for the closure material. Sandover determined 
the mean velocity through the gap using Valentine's Cy as in the first 
procedure. A logarithmic velocity distribution, normal to the curved 
boundary was assumed to obtain the velocity against the advancing face 
as a function of the mean velocity. From this the maximum shear stress 
on the advancing face was determined as another function of the mean 
velocity and was equated to a critical tractive stress for the material 
on the sloping boundary (Carter, 1953). In the analysis, coefficients 
were introduced and were determined empirically as in the first pro- 
cedure. Curves for evaluating the coefficients have been presented. 


The governing equation for stability formulated in this case also needs 
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to be balanced byatrial and error procedure and takes the form 


1/3 at 
Ghobeale) aczuke logs (=) (2.16) 


d D 
in which Kis a constant for one size of closure material and C is an 
empirical coefficient. 
Pariset and Hausser (1959) studied the stability for both 
Straight and angled tip closure. Without elaborating the basis, they 
have developed an empirical relationship for stable closure material 
size (refer Fig. 1.2c) for straight closure as 


é 3/2 


H 
D_ = 9.80 (44) (2.17) 


D8 H 


in which Hy = height of the dam. 

In the final stages of closure, the flow through the opening 
is minimal, when an alternate diversion is functioning. Then both AH 
and Hp approximate to H, which according to equation 2.17 results in 


an extremely uneconomical size of closure material given by 
D = 0.8H (2.18) 


In this procedure, determination of the Stable size needs in 
the first instance evaluation of the head loss at the contraction and the 
upstream backwater depth; both of which can only be evaluated by referring 
to other empirical coefficients. Hence this analysis seems to have 


limited applicability. 
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2.1.4 End-Dump Closures on Alluvial Channels 


No systematic investigation of closures on alluvial channels 
appear to have been reported. Because considerable bed scour occurs at 
high contraction ratios in this type of closure, Izbash (1970) recommends 
rivers with rocky beds be closed preferably by end-dumping. Scanty 
qualitative information is available for end-dump closure on the alluvial 
Volga river (Izbash and Lebedev, 1961). Kuul (1970), Ovchinnikov (1970) 
have reported that the flow demanded an increase of stable size with 
increasing head-differential across the dam. Relevant analysis for 
determining the stable size has not been reported. 

A protective blanket on the alluvial bed is generally provided 
prior to end-dump closure. In such cases if the bed scour is totally 
prevented, analysis for rigid bed channels would be applicable. 

Sandover (1969) has obtained the discharge coefficient of 
end-dump constrictions on a coarse sand bed in a 2 ft. wide flume. These 
coefficients were consistently lower than the corresponding coefficients 
of Kindsvater and Carter (1953) obtained for similar sloping spill-through 
constrictions on rigid beds. Sandover has not advanced any reason for 
the lower coefficients, except rather inconclusively attributing it to 
closure material roughness, although it would be more reasonable to 
expect that the self-formed faces of the closure dam would tend to con- 
tract the flow as a venturi and thereby increase the discharge coefficient. 
In any case, a reasonably accurate determination of either the mean 
velocity or the velocity distribution against the end-dump face is not 


possible with the present knowledge. Sandover (1969) however states 
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without ample justification: “evidently the downstream scour will not 
affect the stability of the dam, unless it is sited on a very fine 

muddy ooze, which might move away from under the dam". Presumably based on 
this, he makes a statement (1971): "this scouring, however, did not appear 
to have any effect on the stability of the dam and consequently this 


work should be applicable to dams founded on loose erodible bed material". 


2.2 Bed Scour Due to End-Dump Constrictions 


2.c.| Mechanics of Scour 
A channel constriction generally results in two kinds of scour. 
(a) Obstruction Scour or Local Scour 
Smal] obstructions such as piers, abutments, or short spur 
dikes affect channel flow only in their vicinity. A vortex system forms 
near the obstacle and the increased bottom velocities associated with 
the vortices increase drag and lift on bed particles, which in turn 
causes local erosion (Neill, 1964, Vinje, 1967). Observations on abut- 
ments and piers by Laursen (1953), on vertical faced and other abut- 
ments by Liu, Chang and Skinner (1961), and on a vertical faced spur 
dike by Garde, Subramanya and Nambudripad (1961), show that deepest 
scour general ly Becurs around the upstream face of the obstacle. 
(b) Contraction Scour or General Scour 
A gradual contraction of a channel over sufficient length est- 
ablishes a new regime of flow and causes a general degradation of the 
bed in the contracted reach, which is termed general scour. The depth of 


scour is then determined from consideration of the increased bed shear 
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23 
in the contracted reach. (Straub, 1936, Laursen, 1960 and 1963, Komura, 
1966). 

(c) Scour at End-Dump Closure Dams 

The scour due to end-dump closure is essentially a combination 
of local scour and general scour. Typical scour patterns obtained exp- 
erimentally in this study and shown in Fig. 2-2 attest to this. 

At low contraction ratios, each wing of an end-dump dam be- 
haves like a short spur-dike built out from one bank or like a projecting 
bridge abutment and creates local scour around its face. Because the 
geometry of the self-formed face of an end-dump dam is different from 
the rigid face of an abutment or a spur-dike, there is an essential 
difference in the scour patterns. Observations on abutments by Liu et 
al (1961), and on groins by Tison (1962), show that maximum scour around 
a vertical-face obstruction is deeper than that around a sloping-face 
or spill-through obstruction. These observations indicate the shape 
of the obstruction to be an important factor in local scour. The stream- 
lined shape of the end-dump face in contrast to a blunt body produces 
a relatively weak wake and contracts the flow gradually in a venturi- 
shape from the dump line to the vena-contracta. This causes the deepest 
scour to move downstream of the dump line to the vena-contracta. The 
scour hole corresponds more or less to the inverted frustrum of an elliptical 
cone. (Fig. 2.2). 

At higher degrees of contraction, the scour holes from the 
opposing sides overlap and a general pot hole scour pattern is observed 
over the entire gap as in a long-contraction. The deepest scour is also 


located downstream at the vena-contracta section. 
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(c) SECTION ALONG 1-1 


FIG. 2.2 SCOUR PATTERNS AT AN END-DUMP CONSTRICTION 
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Location of the maximum scour at the vena-contracta section 
has been observed in the laboratory by Sandover (1970) for end-dump 
constrictions. Deepest scour was also observed downstream of the con- 
striction for spill-through abutments in the studies of Liu, Chang and 
Skinner (1961) and Tutt (1972). 

Although much analytical and experimental work has been done 
on scour, knowledge is still far from complete. The kinematics of 
flow, particularly vortex formation close to obstacles which produce 
local scour, has received little study. Design recommendations have 
been drawn largely from empirical data for specific constrictions. 
Useful summaries on scour at constrictions have been prepared by an 
ASCE Task Committee on Preparation of Sedimentation Manual (1962), by 
Neill (1964) and by the Highway Research Board (1970). 


2.2.2 Maximum Scour Depth 


A distinction is usually made between the "maximum scour 

depth" generally referred to in the context of the clear-water flow 

and the “equilibrium scour depth" generally referred to in the context 
of sediment-transporting flow. Komura (1966) describes these as static 
and dynamic equilibrium scour depth respectively. "Limiting scour depth" 
is also used to denote either of the terms. In this study the term 
"maximum scour depth" is used to denote the maximum depth below the 
original bed level whereas (d.th) represents the maximum scoured depth 


below the original water surface. 


as im" 
: a 


fiohtase sh ait di tity de : oft 
qawh-bns yor iover) VOVOBNES io iapeeel? alent 9 
“nod ans To niserid aniwot bayiedde aS Shia 
bing pend: ats io esthutensae a etapa Re babi beak 
(S0er) douT Hin 
Soh Mead’ Si AoW 16: J namheggs Gna (aakgutais NOt join 
10 2atdenantd oT .stetomos mort vs? diy 2%, spbet wand y 
anuberg dorkw esfuszecn’ oF ‘seala notrambat obra is 
aver enorishisnmossy ipread.. .VouUFY sitant baybegad’ on 
Jnotiorvdanos 15! mVE 7302 %07 6750 foo tortqina mart hagres ay 
5 vd bSTsqStd Need svér 2iprtotnsenos 36 boas : i 
vd ,(S8G!) leuneM norzernantbse to ROrRGaangny, vo deena 
(OSGI) Svsed firoacant va a wi bate: wewne 


ie ts Pie a 
a AM 5 ma ; * 


vudda muntxea? Srtd a tine be ote ie 2 ‘dottont get A 

wort 1s2eweyeat Sendo iibanos, on at of erst upe 

sKadNOD SAS nt Od) baretey el énanap head supa? mieten 

sitste. a6) o2old' gadbisesb | (aaety sun}. Mlk pore 

"addab Wd eee int” ulaiiaegasy Ataab “wage sid tp 9 

mis3 ott ybuse ahha: at eet oil? ‘o wart hs. arpa 63 be 

an woled Agqsb oer sity! stig: oy bee at ons 
dttish. ba wos mint ny sit! 2dndeaiday ie) ‘280% aor, 


aishiid 8 


26 
Assessment of local scour as a consequence of an increase 
in the discharge intensity, has probably been the earliest approach. 
Lacey's studies on scour in alluvial rivers in regime, led to the 
formulation of a relationship by Khosla (1936) between de, the scour 
depth below the original bed and q, the discharge intensity per foot 


width of the main channel, as 
h +d. = 0.90 (estes) (2.19) 


in which f, is Lacey's silt factor, which is a function of grain size 
of the bed sediment. The application of this formula to the calcu- 
lation of scour at channel modifications such as bridge piers, bends 
or spur dikes, requires the use of a coefficient, which ranges from 
1.0 to 3.5 to account for the flow concentration. 

Blench (1957) reported an empirical correlation between 


(htd_) and q developed from prototype data by Andru (1956) in the form 
(reddy Re ean aageq 8 (2.20) 


in which Fe is the bed factor. This is similar to Lacey's analysis. 


Izzard and Bradley (1957) and Ahmad (1953) suggested relation- 


Ships in the form 


2/3 
(h+d.) = K (qo) (2.21) 
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in which Ke a multiplying factor, is primarily a function of the spur- 
dike configuration and qo is the discharge intensity at the constriction, 
(qo = q/1-m). 

The equations 2.24 through 2.27 indicate that the scour depth 
is primarily a function of the increased discharge intensity, which 
is a consequence of the constriction, but the functional relationship 
is complicated, being dependent upon other factors such as the geometry 
of the constriction, the bed material characteristics, the approach flow 
conditions and so on. 

An analysis of scour in a long contraction by Straub (1934) 
and later by Laursen (1960) and Komura (1966) utilising the continuity 
equations for both discharge and sediment, have been based on the 
concept that the scour in a constricted reach would attain equilibrium 
after the sediment inflow balances the sediment outflow. Thus a dis- 
tinction is made between the clear water scour and scour in sediment 
transporting flow. In the case of clear water scour the sediment 
inflow is nil and under equilibrium condition the boundary shear stress 
in the scoured region is the critical shear stress for the material. 

It is obvious that for the same approach flow conditions, the maximum 
scour depth for the clear-water case and the sediment-transporting 
case would be different from each other. Laursen's (1953) experimental 
observations indicate the maximum scour depth in the latter case 

to be largely independent of the flow intensity and size of the bed 
material but dependent principally upon the flow depth and length of 


the encroaching abutment. 


Laursen (1960) considered a sudden constriction due to a 
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bridge as a fore-shortened long contraction" (Fig. 2.3) and proposed a 
relationship for scour at an encroaching abutment under bed-load move- 


ment condition as 


d d 170 


= 2.75 = Kes (A+ - 1] (2.22) 


>|~ 


in which a = length of the encroaching abutment and r is a ratio of 
depth of scour at an abutment to depth of scour in a fictitious long 
contraction of width a + 2.75 d.. 

For the case of clear water scour at an encroaching abutment 


Laursen (1963) proposed 


-2.7 2 (+0 M2) -1] (2.23) 


>| 


in which To and T. are the boundary shear stress in the normal channel 
and the critical shear stress for the bed material respectively. Laursen 
does not consider the effect of the variation of r as significant. In 
Fig. 2.3. design curves to determine d. from plots of equation 2.22 

are drawn for r = 12 (based on experience of sediment-transporting flow) 
and for different values of To/ Te: On Fig. 2.3 equation 2.23 is plotted 
in dashed line to determine d. for sediment-transporting case. These 
curves are valid so long as the scour holes due to the two opposite 
abutments do not overlap. Laursen reasons that when the clear opening 

1s more than 5.5 do; interference is minimal. For further reduction 


in the opening, interference increases and at maximum interference the 
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scour depth is that due to a long contraction. In such a case the 


depth of scour for the clear water case is given by 


(h+d_) Tee 3/% 6/7 
Seer B 
(=) (2.24) 
C 
and for the sediment-transporting flow 
(h+d, ) B 0.59 to 0.69 
= (>) (2.25) 


h b 


Design curves to determine d. from equations 2.24 and 2.25 are shown 
Pee ta.2.4. 

Considering the complexity of the problem, Garde, Subramanya 
and Nambudripad (1961),and Liu, Chang and Skinner (1961) have analysed 
the problem from a dimensional approach. 

Garde, Subramanya and Nambudripad accepted the form of 
correlation of qo and h + d. (equation 2.21), but considered that 
logical dimensionless terms are inherent in the correlation and made 
dimensional analysis of all the pertinent variables. They studied scour 
around a spur-dike made of a vertical steel plate constricting an 
alluvial channel in the laboratory. Bed materials of median diameter 
0.29 mm, 0.45 mm, 1.00 mm and 2.25 mm were used. They obtained a 
systematic correlation between the dimensionless scour depth (h + d.)/h 
and Froude number of the approach flow F and the opening ratio M for 
one size of the bed material. The variation of the three parameters 


for a sand of 0.29 mm was obtained by them as 
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(h+d_) 
= 4.0 (9) Fe/S (2.26) 


They concluded that in general a relationship exists in the 


form 


(h+d_) 
5 = 9 (q) FP (2.27) 


in which J and p are functions of Cy» which is the drag coefficient of 
the bed sediment = UGA S YoY S) d}/we Py They found that J varied in 

the range 2.75 to 5 and p varied in the range 0.67 to 0.90 for the bed 
sediment ranging from 0.29 to 2.25mm. Their experimental data con- 
vincingly indicated that even for sediment-transporting flow the depth 

of scour is not really independent of the flow intensity or the Froude 
number of the flow, although in most of their runs, prior to the introduction 
of the spur-dike either no movement, or only weak bed movement was occur- 
ring. Their data, however, showed a systematic trend when plotted on 
Fig. 2.3 thereby indicating that Laursen's equation 2.28 and 2.29 

have merit. Neill (1962) argued that Garde et al's assumption that Cy 
adequately accounts for the sediment size is incorrect, because beyond 

a median size of 1.5 mm the fall velocity being proportional to the 
square root of size, Cy would be independent of the size and thus depth 
of scour would be same even for large boulders paving the bed. Neill 
(1962) basing on Blench's (1957) bed factor Fie pointed out that an 
analysis of Garde et al's data indicated that the scoured depth might 


be expected to vary as qo 1/8. 
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Liu, Chang and Skinner have conducted extensive laboratory 
studies on the scour at bridge abutments of different shapes such as 
vertical wall, wing-wall, spill-through types. Their findings basically 
confirm those due to Garde et al, that the Froude number of approach 
flow F, as well as the opening ratio M are the most essential non- 
dimensional terms influencing scour. 

For the case of sediment-transporting flow they obtained 


empirical equations for vertical wall models in the form 


(h+d_) F Oden 1/3 
+ =ylUhehs ati 3ili) F (2.28a) 
maximum 
and 
(h+d_) 4 1/3 
| a =2.15 (=) F (2.28b) 


equi librium 


For the clear water scour the equation given was 


d 
| Hi 2 = 12.5 (=) (2.29) 


F 
M 
limiting 
The scour depth for wing-wall and spill-through type abut- 
ments were lower than those for the vertical wall models for the same 
values of F and M. 
Because only two sizes of sand of 0.56 mm and 0.65 mm were 


used in their experiments, no quantitative analysis of the effect of 


the sediment size has been done by Liu, Chang and Skinner. 
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2.3 Maximum Backwater at End-Dump Constrictions 
2.3.1 Constrictions on Rigid Bed Channels 


Lane (1920) is one of the earliest investigators to study 
open channel flow through sharp-edged vertical constrictions. His 
analysis was mainly based upon formulas by D'Aubuisson and Weisbach 
(Yarnell, 1934), where piers formed the obstruction to flow. He did 
indicate that a correlation may exist between the backwater ratio and 
coefficient of discharge, but no definite correlation was given. 

Kindsvater and Carter (1953) and Tracy and Carter (1953) 
did extensive investigations of open channel flow through abutments 
of various shapes and geometry, using the approach of dimensional 
analysis. Their procedure of evaluating the maximum backwater can be 
Summarized as follows. 

Referring to Fig. 1.2b and applying the energy and continuity 


equation, a discharge equation is obtained as 


Q = Cy prgL2a(atagy*/20rhe I) (2.30) 
in which Q = discharge in cfs 
Cy = Kindsvater's discharge coefficient 
b = mean width of the contracted opening 
hy = flow depth at section 33 
y.2 
Oh, a = weighted velocity head at section I-l, where the mean 
velocity is Vy and Oh, is the energy coefficient 
he = head loss in friction between sections I-] and 3-3 
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By the aid of dimensional analysis Cy is found to be primarily 
dependent upon Fs the Froude number of the flow at section 3-3 
(F, = Q/bhatgh,]'/*), m, the contraction ratio and L/b, the ratio of the 
length of the contraction in the flow direction to the opening width. 
Base curves to evaluate Cy correlated to m and L/b, for the case of 
vertical faced constrictions with square-edged abutments placed sym- 
metrically at right angles to flow having Fy = 0.50, have been pre- 
sented. Curves for adjusting the discharge coefficient for other values 
of F3 and other shape and angularity of constrictions have also been 
presented. The objective of the investigation was to predict the dis- 
charge through the constriction on natural rivers where hy and AH are 
known from field observations. 

For similar vertical faced constrictions with square-edged 
abutments, Tracy and Carter have extended this study and have presented 
base curves for a dimensionless backwater rise h*/AH. By dimensional 
analysis h*/AH was empirically correlated with the contraction ratio m 
and the Manning's roughness factor n. Adjustment curves for h*/ AH 
were also presented for other abutment shapes. 

This method cannot be used to estimate directly the maximum 
backwater h* from h*/AH relationships because AH, which is dependent 
upon Q, hy and Cys has to be evaluated in the first instance by a trial 
and error procedure. 

Liu, Bradley and Plate (1957) recognizing the inherent diffi- 
culty in Tracy and Carter's procedure conducted studies for vertical- 


board, wing wall and sloping spill-through abutments. Analysis by 
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36 
energy, momentum and continuity principles led to an equation for H, the 
depth of maximum backwater section, given by 


3 
H 
(th 


> Fook tect vl) (2.31) 
in which F is the Froude number of the uncontracted normal flow; ® is 
an empirical coefficient; M = opening ratio = I-m. 

The empirical coefficient ® was observed to depend upon the 
model type, F and M and is empirically correlated for a simple vertical 


board model as 


® = 1.33 [1 - owe (2-n- 4) (2.32) 
( 3F 


The combination of equations 2.31 and 2.32 led to 


(Fy - 1 = 4.83 FF Fa - £ (2.5 - ") (2.33) 


M 


The ultimate form of presentation of Liu, Bradley, Plate's 


data is a series of curves of (H/h)> plotted against F2(1/Me - 1). 


Their entire experimental data by dimensional analysis led to the 
inter-relation between essential dimensionless variables of the form 


n* = function (— , 2 , model form) (2.34) 


h ¥gh 


wo 


which the analytical relationship in equation 2.33 confirms. 
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The essential feature of their study is that the maximum back- 
water rise h*, primarily a function of the uncontracted flow Froude 
number and the opening ratio, can be directly determined for different 
forms and shapes of constriction from their empirical correlations. 

Studies similar to Liu, Bradley and Plate have been conducted 
for rectangular constriction plates by Valentine (1958) who also adopted 
a dimensional approach. Valentine assumed weir flow through the con- 


traction to result in a discharge equation of the form 


Cee (2.35) 


in which Cy = discharge coefficient, which Valentine determined empiri- 
cally as a function of F and m and presented correlation curves of Cy 
and F with m as a third parameter. 


Ve for V, Q can be written as 


Substituting F(gh) 
Gren eey hemor rig) othe (2.36) 


From equations 2.35 and 2.36 


af é 


(*) = (g/*/C4) - F (2.37a) 
2/3 
or (F) 5 CaS (F) (2.37b) 


Biery and Delleur (1962) studied backwater rise through arch- 
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bridge constrictions. They postulated that equations 2.31 and 2.37b 
are similar, because if the leading term in equation 2.31 is only con- 
Sidered, both the equations indicate a generalised backwater equation 


of the form 
H y 
= Ff) (2.38) 


From their experimental results,Biery and Delleur presented 


an equation for any type of constriction in the form 


3.39 
(F) = 1+ 0.47 (RY?) (2.39) 
in which M' is the representative opening ratio. Liu, Bradley and 


Plate's data fitted equation 2.39 satisfactorily. 


2.3.2 Constrictions on Alluvial Channels 

Although considerable studies have been made to assess the 
maximum backwater on rigid bed channels, the only reported works on 
alluvial channels are those of Liu, Chang and Skinner (1961) and 
Sandover (1969). 

Liu, et al's extensive study of the backwater rise at con- 
strictions in alluvial channels is an extension of Liu, Bradley, and 
Plate's (1957) because Liu, Chang and Skinner accepted some form of 
correlation between the backwater ratio H/h and F°(1/Me - 1) (equation 
2.31) even for alluvial channels, although energy and momentum equations 


for flow over the scoured bed at the constriction would be appreciably 


aint abdutadid peypaee ey ’ . fe 


— 


(9.3) ~ er a | 


boinidesg walled ons qs te aiiwege neal tant 
nat sat wh woregeatenas: +0 9g wm 


(ee.3) a, wo + +8 ote 
| Adee 
bis. yatberd utd ofget eatrieye iain aa : 
yl nnaeneney 
| ie 


oti eese25° 6S sbem vial a) 


ho adiow bet-roqen! yine. od (2 nosy bad bipty. 0 votswilosd 


bas (feet) eal a ozo oe a 


“(105 38 s2l¥ Vat eWAd 
bné . vSlbavt ~wteyer nc 


sere write 2"{s $e wid 
aN 2 ztenseie feivelts at 


40th oF) nb’ ae a ozuaoad <er 
fidhtsups) ~ SHA yP9 bag: d\H oft ortey 


enokinues misasmom bra viene Medonsls pr evs (TELS 
Midst rsnqde od bivow bideideh es ke | 


; 


- 


39 


different from rigid bed flow. They considered AH, the head differential 
to be a better representative term than h*, the backwater rise and com- 
pared their AH/h with that due to Liu, et al (1957) as shown in Fig. 2.5. 
They made a rather qualitative conclusion that "AH for alluvial channels 
is about sixty per cent of that in a rigid channel" 

Using a 2 ft. wide sand-bed flume, Sandover (1969) studied 
the backwater effects of end-dumped dams formed of materials of three 
nominal sizes (D = 0.05 inch, 0.22 inch and 0.44 inch). Sandover emphasized 
the difficulty of applying Carter and Tracy's method involving the un- 
known term AH. He therefore took the dimensional analysis approach 
also and obtained empirically the primary relationship between the 
non-dimensional variables of the form 

V Dies 


h* ’ b 
~—= function (—,%5.,7., >) (2.40) 
h /gh Sa Wee 3 8 


Sandover observed the terms b/L and b/D to be relatively un- 
important and has presented curves for h*/h versus m with V/vYgh as the 
variable. Sandover's observed backwater coefficient h*/h, when com- 
pared with Kindsvater's curves for similar type of sloping wall con- 
strictions, were found to be consistently lower. Interestingly enough, 
Sandover did neither consider the size of bed material as an essential 
independent variable nor mention the effect of bed scour as pertinent, 
although he stated qualitatively that "deep scour hole just downstream 
of the toe of the dam and raised dunes on either side of this hole 


occurred when the larger materials were used as closure material. 
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2.4 Concluding Remarks 

In this chapter, previous studies relevant to the end-dump 
closure of alluvial channels have been reviewed. These studies pertain 
to the stability of closure material both in transverse-dump and end- 
dump closures, the bed scour, and the backwater rise due to constrictions 
in alluvials channels. 

From review of previous research, it is evident that a systematic 
study of end-dump closure of alluvial channels has not been undertaken 
so far and that further study is needed in order e develop comprehensive 


hydraulic design procedures for such closures. 
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CHAPTER III 
THEORETICAL CONSIDERATIONS 


3.1 Stability of Closure Material 


3.1.1 General 

The present state of knowledge on the mobility of a stone end- 
dumped into the flow relies either on the maximum velocity against the 
stone (equation 2.4) or on the mean velocity and depth over the stone 
(equation 2.13). None of these quantities can be precisely determined 
for flow through an end-dump constriction even for a rigid-bed channel 
as clearly indicated in the foregoing discussions. The inaccuracies 
involved in computation of velocity and depth through the gap on alluvial 
Channels would be even greater. A simple velocity approach, with the 
assumption that it accounts adequately for the entire hydrodynamic ef- 


fect of the flow, suffers from a number of deficiencies. 


3.1.2 Deficiencies of the Previous Design Methods 


Detailed information regarding the flow velocity distribution 
through an end-dump constriction is not available. As such the recom- 
mendation of Izbash (1970), that a 20% increase over mean flow 
velocity occurs against the end-dump face cannot be considered as ade- 
quately representating the effect of the accelerating flow. Assuming 
that the maximum velocity against the boundary can possibly be evalu- 
ated empirically, no evidence appears to justify that velocity versus 


competent size relationships derived from criteria for incipient particle 
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movement on the level surface of a transverse-dump dam, are directly 
applicable for flow around an end-dump dam. 

Raudkivi (1967) suggests that "impulse forces which are the 
consequence of turbulent agitation should be more important than mean 
boundary shear for entrainment of grains, whenever the boundary layer 
is turbulent." On a level bed, the turbulent fluctuations of drag 
force as high as twice the mean, occur due to velocity pulses up to 
40% above the mean. Such fluctuations can cause instabilities at 
below-competent velocities. The turbulent velocity pulses, due to the 
accelerating flow on the extremely rough end-dump face are more intense 
than that on a level bed and cause substantial-increase in surface 
drag and lift forces. Hydrodynamic lift,which is not usually considered, 
can contribute substantially to instability as studies by Einstein and 
El Samni (1949) and Benedict and Christensen (1971) indicate. 

The greatest single objection to a simple velocity criterion 
is that curvilinear flow with high local acceleration results in a 
higher surface drag than rectilinear flow. For flow through a curved 
trapezoidal channel bend, Ippen and Drinker (1962) have found boundary 
Shear stress maxima 2.40 times of that on a straight level channel bed. 
Their conclusions for curvilinear flow can be summarised as 

(i) The boundary shear stress increases in intensity as well as 
in areal extent with conditions of increasing stream curva- 
ture presumably due to effects of local acceleration and 
secondary motion in the flow. Locations of shear maxima 
were generally found associated with the course of the 


filament of highest velocity and with the zones of locally 
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accelerated motion. 
(ii) The magnitude of shear maxima is not predictable from head 
loss (in terms of v*/2g i.e. by correlation with velocity) 


or by theory. 


Because flow negotiating an end-dump dam is similar to flow in 
an acutely curved reach of channel, the advancing boundary would experience 
considerable increase of shear stress over the normal that occurs on 
the uncontracted channel boundary. Awazu's (1967) studies indicate 
that the shear stress on the bed of a rigid channel due to a sym- 
metrically constricted flow increases to four times the normal, at a 
contraction ratio of 0.40 for an approach flow of Froude number around 
0.50. 

Relevant laboratory data or. field observations of velocity 
and size are scarce and restricted. 

Sandover's procedure, based primarily on the mean velocity 
consideration and using a trial and error solution offer an acceptable, 
but rather inconvenient design technique for closure of rigid-bed 
channels only. The limitations of Pariset and Hausser‘s procedure 
(1959) have already been stated. 

An attractive and ideal solution would be to abandon the 
correlation of stable size with a maximum velocity, in favour a pre- 
cisely determined maximum shear stress (ASCE Task Committee on Sediment- 
ation, 1966) for coarse, granular closure material. Unfortunately 
it is not possible at present (1972) to predict analytically either the 


shear stress distribution or its maximum value for the accelerating flow. 
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In this study evaluation of the effect of shear stress on material en- 
trainment is attempted by a semi-empirical analysis, although the design 
procedure recommended herein does not require any determination of 
Shear stress. This is followed by dimensional analysis of all the 


relevant variables. 


3.1.3 Theoretical Analysis by Critical Shear Stress Criteria 


It is assumed that the threshold condition is attained after 
the maximum shear stress on the advancing boundary exceeds a critical 
value for the closure material. A semi-empirical analysis for closure 
of a simplified rectangular alluvial channel can be made as follows. 

The uniform approach flow has a discharge Q, a mean velocity 
V, a uniform depth h and a uniform slope i (Fig. 1.2b). The bed 
material has a nominal diameter, d. From Manning's formula for uni- 


form flow 


_ 2) 347 ml 
PSs R, is (32%) 
in which Ke is the Strickler's coefficient, which for fully developed 
turbulence and hydraulically rough boundary can be considered as 


(Constant/d!/®); Ri» the hydraulic radius can be considered equal to h 


for wide channels. 


Equation 3.1 can therefore be written as 
Be aie : 


in which J, is a constant. 
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46 
The shear stress on the bed for the normal uncontracted flow 
T) 1S given by 


Tur, Seria ¢ (3.3) 


Substituting equation 3.2 in 3.3 


a eetieviee dyed’ ee ae 
fo) nt/3 Le h ; 
in which Jo is a constant. 
At a contraction ratio m, the maximum shear stress t » on 


max 
the end-dump face will appreciably exceed T (Ippen and Drinker, 1962). 


Awazu (1967) proposes a relation between the ratio of the 
increased bed shear stress downstream of the contraction to the normal 
Shear stress, and the contraction ratio m, in the form 


109, ( iam) = 1.40 m = 0.021 (3.5) 


t 
Omax 


A similar exponential relationship is assumed for the maximum 


Shear stress on the-end-dump face, which gives 


A ° . 
) j,mtj 


‘6 


in which Jy and Jo are variable indices primarily dependant upon the 


Froude number of the approach flow F and closure material size D. In 
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actual practice the closure material is rarely a unigranular mass. In 
this study, the nominal size of the closure material which is a mixture, 
is taken as the Dos which is defined by the sieve size that passes 75% 
of the material by weight. U.S. Bureau of Reclamation (Lane and Carlson, 
1953) recommends Dae as the representative size for coarse granular 
material. Sandover (1971) also considered Doe as the representative 
SihZel. 


The exponential increase of iT /t_) in relation to m was 


max 
verified in this study, as reported in Section C.2, "Experimental measurement 
of drag on advancing boundary". From equations 3.4 and 3.6 
7s j.mtj 
a eoms es tamale) reike toneth) (3.7) 
Z h 
From studies by Lane and Carlson (1953) the critical tractive stress 


on a sloping bank T, is given by 
b 


2 fs Ve 
iets eu en cOS Um (1-tan"6/tan"o) (3.8) 


in which uP is the critical tractive stress for cohesionless, granular 
material on a horizontal bed, 6, is the angle of inclination of the 
slope, > the angle of repose of the material: 


t. for material of nominal size > 0.03 ft. is given by Shields, 
1936) 


ae 0.06 (¥p-Yy)D (3.9) 


Combining equation 3.8 and 3.9 
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Zany oe Ue 
T. TH=EROR0G (Yp-Yw): cos6(1-tan"6/tan“¢) . OD (3.10) 


“b 
Recent studies by Christensen (1971) indicate that the stable 
size determined by equation 3.10 is somewhat unsafe, because hydrodynamic 
lift has not been considered for the incipient motion. Christensen has 
proposed a correction factor L. which when incorporated in equation 3.10 
yields the correct form of the "size-critical shear stress relationship 
for a sloping bank" as 


2 2 V/2 
TO OF {0.06(y)-y,,) cosé (1l-tan°6/tan“>) .L.}.D (Seb1) 


“b 

L. has a complex functional relationship with the particle 
size Does the roughness of the slope Kes the angle of repose ¢ and the 
inclination of the slope 6. Christensen has proposed design charts 
correlating L. with Dac: Kos o and 6. When the roughness/particle size 
ratio is small, L. can be high enough to cause instability even on ap- 
parently stable flat slopes. L.. becomes insignificantly small as K./Do¢ 
increases. 


Based on the foregoing considerations equation 3.11 can be 


written as 


Tauie 06 Vibear D (3.12) 


in which J. is a coefficient functionally related to 6 and ¢. However 


J. varies within small limits because 6, the angle of inclination of 


3 
the self-formed dam face is fairly constant and the variation of ¢ for 


Bb 
\ 


fies e -anis-tbad' heat bart bbe imag ae 
Simsnxboubyd S2ueoad vocab aF O08 
aol deznsdeisdd  .ao?tom sna tg ont: oft 0 
of :€ Woldeude ni bartering oon maiW sate mt notes? a " 
qiflenottelan e2erte teeta ai al on ha ait ramen 
: eck hia . 


4.0) te ah ‘® nad\@rnés=£) Geqs “Gewnlanle 


atotoven, and idiw qtdanotisrad: ‘Teaat sone ction ae a 

ait DRS ¢ Scagst To afore alt nya agobe anit a) -zeantieuon aft a P 
riedd nphesb béenneyg zen necabsatads # egohe wis te vy - 

oxhe alatd ‘eq\zzonipuoy ori nee alee fxg aw attee 
-96 no neve Ybit idee 92059 ot fovicms fiptd od neg 7 «Theme ato 


sc \,! 26 Trem UCineoi thaphent! asta 2 vouquie aay arden. 


ad feo (Te wotteuas iat Stab ano eat ogarat hie 


— 


a n 


\ 


(SF.€) Q. os e 


fovewdH 26 bas 6 af cate ftaant0r Sogter een 6 Ahi a 
fo nottenttont Yo sfgne sit 0 seumaadh apie Sfame. abittw zotaey gb. 
NOt 70 nottobmny edt Ses suistena Nhe et east mab bonret-tes: et 


49 


material coarser than 0.03 ft. is not appreciable (for example, for a 
typical slope of 2 horizontal to 1 vertical, the term cos6(1-tan@6/tan“4) 
varies from 0.63 to 0.774 for ¢ ranging from 35° to 45°). Considering 
the complex nature of variation of Lys equation 3.12 can be further 


simplified to 
me) 635513) 


in which Jy is a variable coefficient. 
Equating 3.7 to 3.10 for the threshold condition 


1/3 j ym j- 
=e ee (3.14) 


ec 2,d 
J,D = Jo V (7) 

Dividing both sides of equation 3.14 by gh for non-dimension- 
ality and writing F* for V/gh 


1/3 dame j 
ee a re eae (3.15) 


The functional relationship of the terms in equation 3.15 
can only be evaluated by thorough experimental investigation. Assuming 
the variation of Ja to be in a small range, a semi-logarithmic plot 
of (D/h) (h/d) '/3 versus m should be linear for one Froude number of 
approach flow. 

Equation 3.15 demonstrates that the mobility of the closure 
material on the end-dump face is a multidimensional problem. It was 


therefore considered pertinent to apply dimensional analysis to all the 
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variables involving flow, fluid, closure material, and channel bed 
material. At the outset of the experimental programme, in a series 

of preliminary runs, the velocity distribution and the maximum velocity 
occurring against the end-dump face were determined to study the ef- 
fectiveness of the critical velocity criteria. In these runs a dis- 
tinctive nature of the collective movement of the material, when pushed 


into the flow, was also observed. 


3.1.4 Concept of Closure Material Efficiency 


Before attempting dimensional analysis it was considered neces- 
sary to define the stability of closure material not merely qualitatively 
(for example by "first evidence of erosion" by Sandover (1971) or "be- 
ginning of entrainment of material" by Pariset and Hausser (1959) but 
rather quantitatively by introducing a numerical measure termed "effic- 
iency of the closure material". The preliminary runs showed that at 
the threshold condition for a particular closure material the smaller 
sizes in the mixture become entrained rather selectively. Continued 
dumping of the same closure material (to constrict the channel further) 
resulted in the larger sizes becoming entrained at an increased rate 
and carried downstream. The advancing boundary slope flattened and pro- 
gressed in a tongue-shape along the flow. However, some gain in con- 
traction ratio over that corresponding to the threshold condition was 
attained. In an actual closure operation, where a mixture of quarry- 
run rock is pushed, a similar pattern of movement would probably occur. 


Therefore a term, efficiency of closure, is introduced to denote the 
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5] 
percentage by volume of closure material remaining stable within the 
limits of the dump line (i.e. within the useful area), out of the quantity 
duped into the flow at any stage of closure. Figure 3.1 explains this 
with more clarity. 

It is obvious that the efficiency of closure is unity until 
the onset of instability and that it decreases with increase in con- 
traction ratio. The efficiency becomes nil when the closure material 
is completely carried downstream upon dumping and no longer results in 


any advance of the dam face. 


3.1.5 Dimensional Analysis 


The variables entering the problem can be grouped into the 
following categories. 

(1) Geometry of the channel and the end-dump dam: 8B (approach 
channel mean width, b (mean gap width). The geometry of the end-dump face 
is considered constant, because the face of the dam is self-formed. Con- 
striction from one or both banks normal to the flow is considered. 

(2) Flow characteristics: V (approach channel mean velocity), 

h (approach channel mean depth). 


(3) Fluid characteristics: (specific gravity of fluid (water)), 


Pw 
v (kinematic viscosity), g (gravitational acceleration). 

(4) Bed material characteristics: d (median size), Sad (geometric 
standard deviation), P. (specific gravity). 

(5) Closure material characteristics: D (representing size, 


passing 75% by weight), OgD (geometric standard deviation), Pp (specific 
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AT CONTRACTION RATIO m = 
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€, = USEFUL VOLUME OF MATERIAL 
REMAINING IN CLOSURE AREA 


e2e4n 


EFFICIENCY OF CLOSURE, 7 = 


€ = TOTAL VOLUME OF MATERIAL 
DUMPED OVER AREA 12341 


DEFINITION DIAGRAM FOR EFFICIENCY OF CLOSURE MATERIAL 
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53 
gravity (assumed same as bed material)), a(angle of repose in water), 
n (efficiency of closure material). 

As the flow against the coarse closure material will be fully 
rough turbulent, kinematic viscosity v was considered unimportant and 
left out of the analysis. 

Assuming that a particular size of closure material will have 
a particular efficiency of closure corresponding to a particular con- 


stricted width, the following function can be formulated. 
n= ,(B, V; hs os ds ds Sgn? Pp? gps % b) (3.16) 


Assuming repeating variables V, Py and h, application of the 
m theorem yields the following combination of non-dimensional parameters. 
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With the manipulation permissible within the theory of dimensional 
analysis, equation 3.17 can be rewritten as 

ae 
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2h? aie ry Cade °gD? 0.) 


(3.18) 


As this study is limited to experimental closure material of 
specific gravity Py * 2.65 and experimental fluid water, the ratio Pp/ Pw 
will be a constant and is left out of analysis. The terms Ogd? gD? and 


a can be reasonably considered as secondary. As such their effects have 
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not been specifically studied. For these restricted conditions the 
stability equation reduces to 
Pree (3-19) 

The principal non-dimensional terms influencing stability are 

considered to be m, the contraction ratio, V/Ygh, the Froude number of 
approach flow, D/h, the ratio of closure material size to approach flow 
depth and d/h, tne ratio of bed material size to approach flow depth. 
The effect of d/n which is a pertinent non-dimensional term for alluvial 
Channels, iS relatively minor in comparison to m, F and D/h, because the 
term (nya) 43 appears in the stability equation 3-15. B/h is considered 
to have the least influence on stability. 

In previous studies of flow through open channel constrictions, 
by Kindsvater and Carter (1955), Valentine (1958), Liu, Bradley and Plate 
(1957) and Liu, Chang and Skinner (1961), the Froude number of approach 
flow and the contraction ratio were observed to be the most significant 
variables. The ratio D/h in equation 3-19 can be considered as re- 
presentative of the relative stability of the closure material against 
this approach flow. 

Equations 3.15 and 3.19 led to the experimental investigation 
in the first instance, of the relationship between the principal non- 
dimensional terms F, m, D/h and n. The influence of d/h and B/h were 


studied in detail later. 
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55 
3.2 Maximum Scour Depth 
3.2.1 Clear Water Case 

The literature review indicates that the maximum depth of 
either local scour or general scour is usually assessed; 

(1) by relating to the increased discharge intensity (a regime 
type analysis) applying empirical coefficients, or 

(2) from empirical correlations involving relevant non-dimensional 
variables describing the phenomenon; or 

(3) from an analysis of sediment erosion in the constricted region 
by relating to the increased bed shear stress (Laursen 1960, 1963). 
(This approach is strictly valid for a long contraction). 

Each procedure has merit and has proved useful in practical 
applications (Highway Research Board, 1970). The principles of analysis 
involved in one procedure are not fundamentally different from that in 
the other procedures, but are really complimentary (Garde et al, 1961, 
Neill, 1962). 

A regime type analysis requires the correct evaluation of the 
increased discharge intensity, which, however, is not possible for the 
non-uniform flow through a constriction. Use of empirical correlation 
of non-dimensional variables is justified if thorough studies have been 
made for the particular type of constriction, and if these variables 
describe the sediment characteristics completely. Shear stress approach 
can only be considered valid if the flow at the location of the deepest 


scour is two-dimensional. 


Systematic scour studies at an end-dump constriction have 
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56 
not been undertaken so far. It is not possible at present to determine 
the maximum depth of scour by a rigorous analysis of the three-dimensional 
flow at a bluff obstruction. An end-dump dam, however, behaves as a 
reasonably stream-lined obstacle. The contracting flow beyond the dump 
line, visualised by dye streaks in this study, indicated a predominantly 
two-dimensional pattern at all stages of closure. From the scour patterns 
obtained in this study it is obvious that the maximum scour that occurs 
close to the vena-contracta (at A in Fig. 2.2) is more a consequence 
of a two-dimensional flow with a high boundary shear (as in a long 
contraction) than of a vortex (as at a bluff obstruction). 

Based on the foregoing discussions an analysis is made as- 
suming the flow at A to be reasonably two-dimensional. The increased 
bed shear stress at the vena-contracta(Section 33 in Fig. 1.2) is con- 
Sidered primarily responsible for causing the deepest scour. The effect 
of sediment size and density is specifically considered in the analysis. 
For the clear-water case maximum scour depth is attained when the 
boundary shear in the scoured-hole is reduced to the critical tractive 
stress for the bed-sediment. 

It is assumed that one-half of a channel having a normal 
flow width B, depth h and mean velocity V is blocked by an end-dump 
dam of length a. Laursen reasons that the obstructed flow hVa is 
restricted to flow over a width 2.75 d. (Fig. 2.3) in addition to the 
normal flow. Therefore due to contraction of the flow, the discharge 
intensity q3 at A (vena contracta location) would be the maximum to 
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The mean velocity of flow V3 at A, at any time t(t = 0 at 
t 
commencement of scour) is given by 


es = q3/(h+d, - AH, ) (3.20) 


in which d, is the depth of scour and AH, is the head loss at time t. 
Because for subcritical flow AH, —< (htd,), (an assumption equivalent 
to neglecting the kinetic term as suggested by Laursen (1963), equation 


(3.17) can be approximated by 
a = q3/(h + d,) (3.21) 


and after maximum scour is attained, V3 ,» the mean velocity at A is 
e 


given by 
et = q3/(h + d.) (au22} 


For sediment - transporting flow, the equation of resistance can be 


approximated in a power form (Keulegan,1938) 


1/6 


Ree a ee 
V, = 7:66 GO) (3.23) 


in which V, is the shear velocity, R = hydraulic radius ~ depth of flow. 
(Error introduced by this assumption is not significant for the two- 
dimensional flow at A), k. = equivalent roughness (a function of V, 

and der) aval 


From equation 3.23 at any time t 
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and after maximum scour is attained 


—} (3.24) 


in which Vx3 and Vie are the shear velocity at time t and critical 

t 
shear velocity for the bed-sediment respectively. Because yee (h+d.) = 
Voy snede)s division of equation 3-24b by 3.24a yields 


Neg h 2m ¢ 23h (htd.) 73 


( ime 7 |) = ‘thea, y’ (3925) 


in which pee the bed shear stress at time t and cn is the critical 
Shear stress for bea sediment. 

At time t = 0, the depth of scour d, = 0 and the bed shear 
stress is eae which is the maximum attained as a result of the con- 
striction. According to Awazu (1967) with the increase in contraction 
ratio m, T increases exponentially and is related to the boundary 


max 
Shear stress for the uncontracted flow i, in the form 


) jamtj 
SuMaXe eG Re (3,26) 


in which j3 and Ja are coefficients primarily dependent upon F. (Base 


10 instead of e is used for convenience in calculation). 


Substitution of equation 3-26 in 3-25 corresponding to time 


t = 0,. gives 
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(h+d liens 
Bee A ——} (B07) 


Substitution of d for kK. and h for R in equation 3.23 yields 


Wa ee) (3.28) 


Tos assuming Shield's dimensionless critical shear stress 


(Laursen 1963) to be equal to 0.04, is given by 


cae 0.04g Koes) ad (3.29) 


Division of equation 3.29 by 3.28 and substitution of Fo 


for ve/gh yields 


2 EES jomtj (htd.) 7/3 
F h _W_) 10 Sad WF S } 


(s3q) (Gg io, h (3-30) 
Equation 3.30 can be rewritten as 
(htd_) 0:29 9. 0:43 O243Gj2meg,) 
S43 0.85 ;h ae 3° 4 
h 0.695 F (7) Fo 10 foes le) 


According to Chitale (1967) for strict correspondence between 
the power law form in equation 3-23 and the logarithmic velocity dis- 
tribution, the coefficient 7.66 and index 1/6 are not constants, but 


variables in the range 3.2 to 7.66 and 1/12 to 1/3 respectively for 
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60 
different bed regimes. The coefficient 7.66 is eliminated from 
this mathematical development. It is instructive to examine the form 
of the equation 3.31 for the index taking the extreme values of 1/12 


and 1/3. For the index taking a value 1/12, equation 3.31 becomes 


(h+d_) 


0.384 ~. 0.46 0.46(j.mti,) 
= = 0.675 F - 0 an 


U9 7h 
(=) Ko. 


d 


(3E32) 


and for the index taking a value 1/3, equation 3.31 becomes 


On 75° -h 3 Pw 0.38 0.38(j.m+j,) (S233) 
(27) (—— 
Ap, 


Because equations 3.32 and 3.33 are not significantly dif- 
ferent from equation 3.31 and the law of resistance of sediment - 
transporting flow is not definitely established, equation 3.31 is con- 
Sidered appropriate for experimental investigation. 

Considering the simplified assumptions made in this analysis 
it is expected that the indices in equation 3.31 would require some 


modification to comply with the experimental data. 


3.2.2 Sediment-Transporting Case 


The approach channel would transport bed load if cae exceeds 
unity. In reality somemovement commences earlier, around a ae value 
in the range of 0.5 to 0.8, as Laursen's plots (Fig. 2.3 ) show. Several 
factors such as grain-size distribution of bed sediment, the variability 


of the critical tractive force particularly for fine sediments, and flow 
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characteristics are responsible for this. In this case, after intro- 
duction of the constriction, the flow as well the sediment input from 
the upstream from a channel width equal to a, (the length of the dam 
Fig. 2.3 ) would be carried over the scoured area in addition to 
the normal flow and sediment load. In other words, if q. is the sedi- 
ment discharge for the normal flow, then thoy the sediment discharge out 
of the scour nole, after introduction of the constriction is given by 


q,. = 9,(1+2) (3-34) 


3 
in which Z, a coefficient, can reasonably be assumed to be a function 
of F and m. 


Equation 3.34 can therefore be written as 


gee ge f,(F sm) (3.35) 


3 

At A, corresponding to the equilibrium scour condition, the mean velocity 

V, and depth h + do, would just carry 4s, such that sediment output from 
e 

the scour hole just balances the input. According to sediment discharge 

formula of Kalinske (1949) 
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62 
in which Vee is the bed shear velocity in the scour hole at equilibrium. 


e 
From equations 3-36 and 3-37 


ees)” 


si a, (3538) 


At time t = 0, (commencement of scour), V, , the bed shear velocity is 
3 


related to the normal shear velocity by equation 3.26. 


=< 


* jam+j, 1/2 
Poe rie as (3.39) 
* 


Application of equation 3.23 to the commencement and equilibrium 


Scour conditions yield 


3 
y= a (3.40) 


Because q3 = V3.h = V3,(htd.), equation 3.40 can be rewritten as 


= (3241) 


Eliminating V, between equations 3.39 and 3.41 and equating 
3 


to equation 3.38 
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Substituting 3.35 in 3.42 and algebraically manipulating 


(h+d_) 0.43(j amtj 
r—— = 10 4 8 re (Fam (3.43) 


Equation 3.40 indicates that [(h+d.)/h] is a unique function of F and m 
and the functional relationship can only be established by experimental 
investigation. 

Equation 3.43 also substantiates the hypothesis postulated by 
Laursen (1953) that the sediment size is unimportant so far as the 
limiting scour depth in flows carrying appreciable sediment change is 


concerned. 


3.2.3 Dimensional Analysis of Maximum Scour 


The variables influencing maximum scour depth d. can be grouped 


into the following categories: 


(1) Flow Properties: V (approach channel mean velocity), h (ap- 
proach channel mean depth). 

(2) Geometry of the channel and the end-dump dam: B (approach 
Channel mean width), b (mean gap width). Self-formed face of the dam 
is assumed to have a constant geometry. 

(3) Fluid Properties: Pw (specific gravity of water), v (kinematic 
viscosity), g (gravitional acceleration). 

(4) Bed Material Properties: d (median size), Sad (geometric 
standard deviation) , Pe (specific gravity), w (fall velocity to account 


for particle shape). 
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64 
(5) Time of scour t: Closure material forming the dam is con- 
sidered to have no significant effect on scour and is therefore left 
out of the analysis. 
Assuming a particular maximum scour depth is attained at a 
particular time t, corresponding to a particular b, for particular 


approach flow conditions, 


d, = F(t, Bs Vs hs Oys Vs Gs ds Oye Ops ws b) (3.44) 


Assuming V, Pw and h as repeating variables, m-theorem 
yields. 


a 


2 0 
is V V d W Vh 
Vie eee ) (3.45) 


By manipulation of the dimensionless groups equation 3.45 


can be written as 


QP 
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Sd has a definite effect on the depth of scour, as studies by Komura 
(1966) indicate, but because his experimental data showed the depth 
ena as Cee for the sediment-transporting case, and cee 
for the clear-water case, variation in gd over practical limits (1 to 
2) would not be expected to significantly affect d.. Hence in this 


study effect of variation of Oud is not specifically studied. Reynold's 


number of the flow Vh/v is considered to have a secondary influence for 
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the turbulent flow in the low viscosity fluid (water). Variation of 
B/h is similarly considered to have insignificant effect. Although 
the fall velocity has been recognised by several investigators (Rouse, 
1939) to be a characteristic parameter describing the shape and erosion 
characteristics of coarse sediments, it is established (Alger and 
Simons, 1968) that variability of shape can considerably affect the 
fall velocity. Correlations with equivalent spherical diameter, as is 
normally done are not reliable. From these considerations, the terms 
d/h and P./Pw are retained and V/w is left out of the analysis. 

For these restricted conditions equation 3-43 can be re- 


written as 


) (3.47) 


From equation 3.47, the functional relationship for the 


maximum scour depth can be written as 


) (3.48) 


in which t. is a characteristic time at which the scour practically 
ceases. It has been recognised (Neill 1964) that te for sediment - 
transporting flow is finite whereas for strictly clear-water flow, 

the limiting scour is attained asymptotically at infinite time. It 
is however reasonably correct to assume the non-dimensional scoured 


depth (htd_)/h to be functionally related to m, F, d/h and p/P, only 
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66 
for both the clear water and sediment-transporting case, whereas a 
correct functional form of equation 3.47 for the clear water case 
should yield (htd,)/h = ] at t = 0 and (h+d,)/h = (htd.)/h at t =, 
(Govind Rao and Sharma, 1967). For the sediment transporting case the 
relationship should yield (h+d,)/h =] at, t.="0) and (h+d,)/h = 


(htd_)/h atit.= to: 


3.3 Maximum Backwater Rise 
3.3.1 General 

Analysis for rigid bed channels by Liu et al (1957), Valentine 
(1958), and Biery and Delleur (1962) have shown that the generalized 


backwater equation is of the form 
H > (e 
he F(a) (3.49) 


Unlike the rigid bed case, backwater rise in the case of a 
mobile bed is influenced by bed-scour. The bed scour downstream of the 
constriction (Fig. 2.2) is similar to a drop introduced in a sub- 
critical flow. It has been established from energy considerations 
that for subcritical flow through a constriction, a drop in the bed 
would result in a higher tail water or less head loss and less back- 
water than a flat bed. Because scour is related to the flow, fluid and 
sediment characteristics, the head loss or the backwater rise given by 
equation 3.49 would logically involve additional non-dimensional terms 
involving sediment size. The closure material size is also important, 


particularly in the case of dams with wide crest and base widths. The 
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67 
effect of turbulent seepage when employing large irregular quarried 
blocks would also be appreciable on the backwater rise. It is thus 
obvious that any generalized backwater equation for end-dump closures 
involving a number of variables would be deficient and cannot be con- 


Sidered applicable to all situations. 


3.3.2 Dimensional Analysis of Backwater Rise 


The variables considered are the flow, fluid and sediment 
properties dealt with in section 3.1.5and 3.2.3,as well as the depth of 
scour d. and length of the dam L, along the flow. Dimensional analysis 


yields a relationship between the most essential variables in the form 


2 >, shee x) (3.50) 


>0 
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It has been established in equation 3.31 that d/h is a function of F, 
mand d/h. In the limited scope of experimental investigations L/h 
variation will not be significant and hence left out. Sandover's (1970) 
study supports this hypothesis. Similarly B/h and D/h are considered 
of secondary importance. 


Equation 3.50 can therefore be rewritten as 
see qd 
nite Fo(F; ms h ) (3.51) 


From the foregoing discussions, it is apparent that the 


experimental data can be analysed to develop correlations( 1) between 
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(H/h)° and Fo(1/Me - 1} according to the findings of Liu and Bradley 
and Plate (equation 2.19) and (2) between C and F with mas a third 
parameter, as suggested by Valentine (equation 2.23). 

In both these analyses, the effect of d/h is critically examined 
to evaluate its effect as equation 3.51 indicates. 

Because both clear-water flow and sediment-transporting 
flow will have the same essential non-dimensional terms given by 


equation 3.51, separate analysis for each case iS not warranted. 
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CHAPTER IV 


EXPERIMENTAL EQUIPMENT, MATERIALS AND PROCEDURE 


4.1 Equipment 


4.1.1 General 
The equipment used in the testing may be divided into 
(1) the flume and the flow system, and 


(2) the measuring devices. 


4.1.2 Flume and Flow System 


(a) Main Flume 

Experiments were carried out in a 60 ft. long non-recirculating 
flume (referred to as flume), 7.5 ft. wide by 4 ft. deep with a hori- 
zontal bottom resting on the laboratory floor. The side walls were 
built of 8 in. wide hollow-concrete blocks. Of the 4 ft. depth avail- 
able in the flume, the bottom 1.5 ft. was filled with the bed material 
over 50 ft. length. A head tank, an upstream stilling pool between the 
head tank and the flume, and a drop box downstream of the tail gate 
took up an overall length of 10 ft. 

The downstream end of the flume was designated as station 0 
and the upstream end as station 50. Of the 50 ft. length of alluvial 
bed, a 30 ft. stretch between stations 7 and 37, where uniform water 
Surface slope could be attained was used as the test area. The bed 


elevation (top of the alluvial bed) was arbitrarily designated as 
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1.500 ft. The bed fill of 1.5 ft. allowed for full vertical develop- 
ment of the scouring action around the end-dump dams, which were built 
symmetrically from both walls at station 26 for all the runs in series 
A, B, C and D and at station 25 in the other series. (Section 4.3.1). 
Flow in the system to the extent of 3 cfs was provided by 
two centrifugal pumps, which drew supply from an underlying sump. 
The flow was received in a head tank and the discharge was measured 
by calibrating the overflow weir in the head tank. Baffles and wire 
mesh screens were provided in the head tank and in the upstream still- 
ing pool of the flume. The downstream end of the flume was provided 
with a tail gate which could be cranked up or down about its bottom 
(horizontal axis) kept at the bed elevation. This tail gate regulated 
the stage-discharge pattern of the normal flow required in the study. 
The flume had provision for a 2 ft. wide rectangular notch 
in the left wall at station 35 for effecting diversion of the flow 
during closure studies. Provision was made to modify the notch to a 
triangular shape if required. The bottom of the notch was kept at 


the bed elevation. 


(b) Subsidiary Flume 


Studies on stability, scour and back water were carried out 
in the main flume. To supplement the studies on scour, a few check 
runs were carried out with pebbles and coal as bed material (Section 
Ave. Vb) A small, 1 foot wide glass walled flume was used for this 


purpose. 
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4.1.3 Measuring Devices 


The principal recording instruments comprised; 

(1) A point gauge for measuring the water-surface elevation and 
the static bed elevation at any point within the flume. The point 
gauge was equipped with a vernier to measure to the nearest 0.001 ft. 

(2) An electronically sensed propeller meter, for velocity 
measurement. The propeller meter was mounted on a special swivel, so 
that continuous velocity traverses could be taken at any inclination, 
particularly normal to the curved advancing boundary. The propeller 
meter was also equipped with a vernier to measure flow depth to the 
nearest 0.001 ft. The propeller meter was connected to an autograph 
recorder. This permitted measurement of velocity to the nearest 0.01 
fps. 

(3) An electronic bed probe for recording bed elevation under 
water as scour progressed. The probe maintained a constant distance 
above the bed and worked on an optical principle. Light was fed through 
optical fibres and reflected back from the bed into a photo cell, 
which operated a servo-mechanism. The bed probe was also connected 
to an autograph recorder and could be calibrated to read to the nearest 
0.001 ft. 

The point gauge, the propeller meter and the bed probe were 
mounted on a single carriage, which could travel along the flume on 
the rails fixed on the side walls. 

Other recording instruments were two steel tapes, fixed on the 


two flume walls and a third steel tape fixed on the carriage to provide 
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readings for the longitudinal and transverse locations of the principal 
instruments. 

A stationary point gauge was provided in the head tank to 
measure the head on the overflow weir. Another stationary point gauge 
was provided at the tail gate to record the tail water elevation. A 
third point gauge was provided at the diversion notch. 

A view of the main flume looking upstream, (Fig. 4.1) shows 
the carriage with the point gauge, the propeller meter, the bed probe 


and the autograph recorders. 


4.2 Materials 
4.2.1 Bed Material 
(a) Sand 
Three sizes of sand of 1.20 mm, 0.60 mm and 0.25 mm median 
diameter, were used as bed materials to cover the applicability of 
this study to closure of alluvial channels. (U.S. Bureau of Reclamation 
recommends the size passing 75% of the material by weight as the re- 
presentative size for coarse non-cohesive materials above 0.2 inch. 
For sizes up to 5 mm, it was arbitrarily decided to accept the median 
size and above 5 mm, the 75% size as the representative size in this 
study). 
The size distribution curves are shown in Fig. 4.2. The 
characteristics of bed material are presented in Table 4.1. 
(b) Other Bed Materials 
The design recommendations of this study are based on the 


results of tests (runs 1 to 100) carried out in the main flume, with 
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TABLE 4.1 
CHARACTERISTICS OF BED MATERIAL 


a 


a 
Material d59 Ps Sgd 
Median Size Specific Geometric Remarks 
in mm Gravity Standard 
Deviation 

Sand T2c0 2.66 WAS Natural river sand 
Sand 0.60 2.65 52 Artificially sorted 
Sand Ge25 S07, lec, Artificially sorted 
Pebbles D>5=6.60 2.68 J 251 Moderately angular 

shape (used only in 

check runs) 
Coal 1.60 eco Crushed coal (used 


only in check runs) 


Soa the geometric standard deviation is computed as 5 (dge/dgntden/dy 6} 


(ASCE Task Committee on Preparation of Sedimentation Manual, Committee 
on Sedimentation, 1962) 


the three sizes of sand as bed material. In section 3.2 the depth 
of scour has been derived as a complex function involving the size as 
well as the specific gravity of bed material. It was therefore con- 
Sidered essential to test the applicability of the design curves in 
a few check runs to 
(i) coarser bed material and 
(ii) light weight bed material. 


The subsidiary flume was used for these runs. 
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Naturally worn pebble having do. of 6.60 mm was used as 
coarse bed material in the check runs. This material was also used 
as closure material B (section 4.2.2) and its size distribution curve 
is shown on Fig. 4.3. Its specific gravity was close to 2.65. 
Coal (anthracite) having de of 1.60 mm and specific gravity 
of 1.26 was used as the light weight bed material in the check runs. 


Size distribution curve for coal is shown in Rig. 4.2. 


4.2.2 Closure Material 

Three sizes of closure material having Do. of 6.60 mm (0.021 “ftp 
18.00mmm (0.059 ft.) and 25 mm (0.082 ft.) were used. These were de- 
Signated as material B, C and D respectively. Closure material A was 
the natural river sand, of 1.20 mm median size, used as bed material. 
This was used only in four runs, for which the approach flow had a low 
Froude number of 0.10. 

The size distribution curves are shown on Fig. 4.3. Figure 
4.4 shows the materials B, C and D in close view. Their characteristics 
are presented in Table 4.2. 

The angular shape of closure materials used in this study 
were considered representative of the natural quarried rocks used in 
river closures. It is obvious that extreme angularity or extreme round- 
ness would change the material behaviour to a certain extent, which, 
however, may not be considerable, because » varies in the range of 38° 
to 42° for variation from extremely rounded to extremely angular shape 


(U.S. Bureau of Reclamation, Sandover, 1971) in respect of materials 
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MATERIAL B 


MATERIAL D 


MATERIAL C 


Do¢ = 25.00 mm 


Dorp = 18.00 mm 


Doe = 6.60 mm 


FIG. 4.4 CLOSE VIEW OF THE CLOSURE MATERIALS 
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TABLE 4.2 
CHARACTERISTICS OF CLOSURE MATERIAL 


Bo Pp gD 9 
Material rai Specific Geometric Shape Angle of 
Gravity Standard Repose 
Deviation in degrees 
A. Sand deol .20 2.66 1.79 Slightly 28-30 
angular 
B. Naturally 
worn 6.60 2.68 eft hema tae 31 
pebble | 9 
C. Crushed Moderately 
Chippings 18.00 2.64 53 angular 34 
D. Pebbles 
and : Slightly 
chippings 25.00 2.66 estZ angular 38 
crushed 


in the size range of 2 in. to 2 ft. Sandover (1971) also observed that 
the angle of repose in water was consistently less than 9 by 10°. 

These observations do not imply that no advantage would be 
gained by the use of extremely jagged rocks, because such rocks by 
virtue of interlocking under the action of gravity and hydrodynamic 
effect of flow, tend to gain added stability. 

This study attempts to make design recommendations in respect 


of closure materials of standard shape and specific gravity around 2.65. 
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4.3 Procedure 
4.3.1 Arrangement of Test Series 
In order to gain an insight into the role played by velocity 
in entrainment of closure material, the experimental programme at the 
outset comprised a series of preliminary runs in which the flow pattern, 
velocity distribution against the closure material and though the 
constriction were thoroughly studied. 
The subsequent sixteen series of runs were conducted to 
yield data for empirical correlation of essential non-dimensional terms 
in respect of stability (section 3.1.3 to 3.1.5 ), scour (section 3.2.1 and 
3.2.2) and backwater rise (section 3.3.2). These runs, covered the 
bulk of the experimental programme. The last four series of runs were 
check runs, meant to 
(1) assess the effect of rate of dumping on stability, 
(2) assess the applicability of design curves for stability to 
irregular channel cross-sections, 
(3) assess the applicability of the design curves on scour to 
(a) coarse bed material (small values of the ratio of depth of flow 
to bed material size, h/d) and (b) light weight bed material. 
It has been shown in Chapter III, that dimensional analysis 
of independent variables describing stability, scour or backwater, 
indicated the Froude number of approach flow to be a significant non- 


dimensional parameter. 
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(a) Froude Number of Approach Flow 


It was decided to conduct stability studies at three repre- 
sentative Froude numbers of 0.50, 0.29 and 0.10, on the consideration 
that in practical cases where closure is attempted, the Froude number 
rarely exceeds 0.30 and is mostly in the range 0.05 to 0.25. Table 
4.3 summarizes the Froude number of normal flow in some rivers during 


actual closure operations. 


TABLE 4.3 
APPROACH FLOW DATA OF SOME RIVERS AT CLOSURE 


B h Q F 
Name of Width Nominal Discharge in Nominal Froude 
River in feet Depth of Cubic ft/sec. Number of 
Flow Flow 

Missouri at a 
Fort. Randal 1,000 12 28,600 Os t2Z 
Missouri at 
Ganeanamessi te 650 24 30,000 0.07 
SOURCES 520 «50 108 ,000 0.103 
Dalles Dam 
River in USSR 158 
(Sandover 1971) 984 27.50 127,000 0 
Mangla Dam 470 18.00 7 ,000 0.07 
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An outline of the experimental programme as split to different 
series is presented in Table 4.4. It can be seen from Table 4.4 and 
Table D.1 that each series (except the check runs) is a convenient 
group of runs corresponding to one Froude number of approach flow and 
one size of bed material. Thus concurrent variation of Froude number 
and bed material could be studied. 

In general, corresponding to the Froude numbers of 0.50, 0.29 
and 0.10, the uncontracted flow had normal depths of 0.18 ft., 0.26 ft. 
and 0.38 ft. respectively. These depths were chosen in order to limit 


the maximum scour depth to 1.50 ft. at the highest contraction ratio. 


4.3.2 Data Taken 
Experimental procedure involved gathering data for varying 
approach flow conditions, i.e. at the three Froude numbers of 0.50, 
0.29 and 0.10, in respect of 
(1) flow pattern and velocity distribution through constriction, 
(2) stability of each size of closure material by obtaining 
closure efficiency - contraction ratio relationships, 
(3) maximum depth, pattern and time history of scour at varying 
contraction ratios, 
(4) maximum backwater rise and its location. 
The general test procedure in each series was as follows. 
The bed was screeded level in the test reach and an equili- 


brium steady flow condition was established by regulating the tail gate to 
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result in the desired depth of flow corresponding to the discharge and 
the Froude number. A uniform water surface slope in the test area 
was attained within a reasonable period of time, but the bed slope 
did not become parallel to the water surface slope, nor could this 
be achieved by shaping the alluvial bed prior to the run. As shown 
in section A.3, the theoretical bed slope for uniform flow was too 
small (0.002) to be maintained precisely in a small reach of 30 ft. 
Thus the normal flow referred to in this study, has a uniform 
water surface slope over a horizontal bed. The slightly accelerating 
flow did not result in any significant increase of velocity to alter 
the Froude number of the flow over the test area, upon which the 
analysis is primarily based. Uniform water surface slope over a 
horizontal bed was also the case in the investigations of Chang (1939) 
and Liu, Chang and Skinner (1961). It is, however, true that open 
channel flow through a constriction being predominantly a local 
phenomenon, the measured data (dependent variables) would not in- 
volve appreciable error with the slight difference in slope of 
the bed and water surface. Simons, Richardson and Nordin, Jr. (1965) 
have demonstrated that normal flow over a sedimentary bed can only 
be considered as being in equilibrium and not strictly uniform. The 
characteristic features of normal flow are detailed in Appendix A. 

A series generally consisted of three sets of runs corres- 
ponding to the three sizes of closure material. The procedure for 
a set of runs was to keep the velocity and depth of the approach flow 


unchanged, thereby maintaining one Froude number of approach flow for 
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all runs, but increasing the contraction ratio from run to run. The 
closure material was dumped from both sides of the flume in equal 
quantities until the dam came up above water level and advanced into 
the flow. 

The rate of dumping and the crest width of the dam approxi- 
mately simulated the following typical prototype conditions: 

(i) 20 cubic yard dump trucks assumed dumping a load every 
five minutes, 
(ii) crest width 40 ft., 
(iii) linear scale ratio 1:100, time scale ratio 1:10. 

At an arbitrarily chosen contraction ratio, the dumping was 
stopped. The water surface levels, live stream boundary, and velocity 
profiles normal to the bed and advancing boundary slope were then re- 
corded, starting at the section of maximum backwater upstream and 
ending where the flow widens to the parent channel] width (section 11 
to section 44 in Fig. 1.2a. The readings were taken fairly close in 
plan around the advancing boundary, and at one-foot intervals else- 
where. The bed profiles were also recorded by the probe at suitable 
intervals to complete the bed configuration, upon which the water 
surface configuration and velocity pattern could be superimposed 
to provide a detailed flow description. 

The configuration of the crest above water and the toe 
location of the end-dump dam were carefully recorded, for computing 


the contraction ratio. These data were used to determine the volume 
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of material moved beyond the useful closure area required for computation 
of the closure efficiency (Fig. 3.1) corresponding to the particular 
contraction ratio. 

By injecting dye from upstream and also using confetti, the 
flow pattern and eddy zones could be visualized. Photographs of the 
flow pattern, the closure material movement, the scour pattern and other 
pertinent features were taken during the tests. 

The flow was continued until nearly equilibrium scour conditions 
were attained downstream. The scour holes were contoured in the dry for 
taking photographs and the scoured bed elevations were recorded. 

The complete test data are presented in a tabular form in 


Appendix D. 
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CHAPTER V 
FLOW CHARACTERISTICS AND MATERIAL STABILITY: 
DATA ANALYSIS 


5.1 Introduction 

Based on the analysis of experimental data of this investi- 
gation, the characteristics of flow through an end-dump constriction, 
and the related problem of closure material stability are studied in 
this chapter. 

The theoretical analysis in Chapter III indicate the con- 
traction ratio, m, to be an important parameter in the hydraulic 
design. Due consideration was therefore given in defining and evalu- 


ating m, particularly for the closure of alluvial channels. 


5.1.1 Computation of Contraction Ratio 


In the previous studies on flow through constrictions, the 
contraction ratio, m, expressed in terms of width, (m=1-B/b) has been 
accepted as the representative measure of channel contraction. For 
closure of rigid-bed channels the representative b is generally evalu- 
ated (Izbash et al, 1970) as 5 (b,+b,)s by and by, representing the 
water surface width and bed width respectively at the dump line 
(Fig. 1.2c). To account for the bed scour in the case of closure of 
alluvial channels, the representative b should be evaluated as 


; b,. representing the width at the scoured bed level. How- 


(b, +b, <) 3 bs 


ro| — 


nt Yoo iy wall ag ames recut "6 sd“ad Mm orton na ei 
A 
vo boa prtnitsb At now tg oF clin ae nol terayen os aad at ner 


insta Tetwwets Yo sized ona 1 Not eit of 


riaiytenoo oud wort 10 gotbuse: iuo)vora att at 
d cad (dN@=Pem) .ddblin Ye aire ait ezeorgne vt cotter We seal 
notsosrtnos Ionnsda Ye swessn ovtssinazorg sone 
noo at d outdsmns2e790% ony tana bed Mota re 
ont tieee 190" atta bk, As om (over ls: yo 28 st) 
b ott $6 eiowtsoeqer ont bas btw sostw 
to Gweels to saan ots — Mooe bee ate : 
ietsufeve ae od Divo 

Pagal 


me hae YY 
- mi) Hy - 7 _ ue 7 
mn . _| _ 
ay a ine pre ; 
. Pol ; 
[ be 2 a 
; ‘a ay ; 
a Or 
: ¥ ; a) Ne h 
‘ iad 
an 
an as hd 
i een _ sy 1G | 
y RATA “1 m : 
: : ia Ui ee : 
. ALATA i ROTTOIRSTOARAH) WOR WOH 
ee bl Sead a lve 
: a i f" ru =) 
peYIAMA ATA 
Ay at ag rn f s Y 
™~ ae ; ; ~, Sian 
| | not?aubo"s 
1 | goes 
isinonitagne Po zie! ans sad tw. ‘bee : 


t 


ie Hosts me 79 ate eras 7 ~~ 
dete falystom one of 0 melding a 


a eee vedas a 


=~ 


Ith yesq BA nt ehgytens. feats round: an 


mu’ ¥ 
J 


ghia olson neo 70 ot nova 2 


al 


Ne F 
say fi 
uh 


a 


ee 


88 


0 0.2 0.4 0.6 0.8 1.0 


CONTRACTION RATIO m 


FIG. 5.1 PLOT OF ACTUAL CONTRACTION RATIO 
m. VERSUS CONTRACTION RATIO m 


ACTUAL CONTRACTION RATIO mg 


ever, Dis cannot be determined precisely for the variable scour con- 
ditions. In their study of bridge constrictions in alluvial channels, 
Lin et al (1961) have stated that "b = width of opening in constriction", 
but did not clarify as to how b was measured in their study. 

For reasons stated above, it was considered practical to 
evaluate b in this study as 5 (b,+b,), irrespective of the bed scour 
at the dump line. The contraction ratio m, is computed as (1-[b, +b, ]/2B). 
It is obvious that the actual contraction ratio ma» (m, =]- [b, +b, .1/2B) 
for a scoured bed is slightly more than the computed value because Dis 
is less than bb: 

From the measured values of by b, and bh. in all the runs 
in series B, C and D of this study a plot of m, versus m_ is ae 


a 
in Fig. 5.1. In spite of the considerable bed scour in the dump area, 
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89 
particularly at higher values of m, the scatter is slight and the cor- 
relation is good. This illustrates that the consistent use of m, 
as the representative contraction ratio, in preference to ma does not 


involve any significant error. 


9.2 Flow Pattern Through End-Dump Constriction 
5.2.1 General 


The longitudinal water surface variation along the centre 
line of the flume at a few typical contraction ratios are shown in Fig. 
5.2. The water surface topography for a contraction ratio of 0.56, is 
aso shown injFig. 5.2. -Rapid drop of water surface contours around 
the dump line in contrast to a flatter drop along the centre line 
indicates the flow to be accelerating more rapidly close to the advan- 
cing boundary. The patterns as visualised by dye streaks are shown 
UE ANC Pe sien c 

It was generally noted that the flow through an open channel 
constriction formed by end-dump dams on a mobile bed is quite similar 
to that on a rigid bed, as observed in the extensive experimental in- 
vestigations by Kindsvater and Carter (1955) for various constriction 
geometries and by Biery and Delleur (1962) for semi-circular arch 
bridge constrictions. 

The progressive closure by end-dumping resulted in a pro- 
gressive rise of the backwater upstream of the dump line. The back- 
water level, which attained a maximum at the instant the dumping 


stopped, was observed to fall gradually, with progression of bed scour 
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(a) At a low contraction ratio of 0.59 
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(b) At a high contraction ratio of 0.82 


FIG. 5.3 FLOW PATTERNS THROUGH END-DUMP CONSTRICTION 
AS VISUALISED BY CONFETTI 
(M in photographs represent contraction ratio) 
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92 
at the contraction, until both attained equilibrium. 

The slope of the dam advancing into the flow was self-formed, 
as an interaction between the flow and gravity. Until the onset of 
instability, the slope more or less conforms to the angle of repose 
in water for the closure material, but beyond the threshold condition 
the slope flattens and the tongue at the end progressively extends 
along the flow. The general configuration of the advancing face and 
scour hole are shown in Fig. 5.4 for two typical contraction ratios 
of 0.59 and 0.77. At the contraction ratio of 0.77, it can be seen 
that the closure material has crossed the threshold condition and the 
efficiency of closure was 0.50. 

The upstream face of the advancing boundary more or less con- 
forms to the contracting live stream. Instability of the dumped ma- 
terial was observed to occur predominantly over the upstream face 
(area pqrp shown cross-hatched in Fig. 1.2a). The balance area down- 
stream on the advancing boundary is not subjected to the erosive 
action of the flow. This was also observed in the laboratory by 
Pariset and Hausser (1959) and Sandover (1969) and in the field in 
the actual closure of Dalles dam (1965). However, against Pariset 
and Hausser's observation that only the upstream 15% of the total width 
of the progressing face is subjected to severe erosion by flow, the 


present experiments indicated the proportion to be at least 60%. 


5.2.2 Contraction and Expansion of the Live Stream Boundary 


The flow contracts from section 1-] through section3-3(Fig. 5.5) 
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(b) At a high contraction ratio m=0.77, F=0.50, h=0.18 ft., d.=0.75 ft. 


FIG. 5.4 GENERAL VIEW OF THE END-DUMP DAM AND SCOUR HOLES 
(M in photographs represent the contraction 
ratio, bed contour interval 0.10 ft.) 
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and accelerates as it negotiates the end-dump boundary. In contrast 
to a relatively shorter length of aren se ee the live stream ex- 
pands gradually over a much longer length. The local angle of con- 
traction Vie at the dump line was observed to be larger than the local 
angle of expansion Pog Measured slightly downstream of the venateontractn 

The flow width at vena-contracta Dos was measured by intro- 
ducing dye from the separation point along the advancing boundary. 
Because the boundary of the jet was not clearly defined, the data 
can not be considered as precise. The length of contraction L53 and 
the length of expansion L34 were also measured. The complete experi- 
mental data are presented in Table D.1. The coefficient of contraction 
C. was computed as the ratio b/d. 

In a rectangular constriction, the contracting stream main- 
tains a more or less vertical edge. For the end-dump constriction, 
the contracting stream edge is more or less vertical on the upstream 
and downstream of the dam and hugs to the sloping end-dump face in- 
between. It is therefore likely that the geometry of the constriction 
has a decisive influence on the coefficient of contraction. Other 
factors upon which C. may depend are the roughness of the closure 
and bed material and the B/h ratio of the uncontracted flow. 

A logical relationship of principal non-dimensional terms 


denoting flow contraction can be assumed as 


oa f, (m,F tang) (oe) 


tenttnon AT .evebRuOd. a ub no 8 
“xe mbo"1S2: VAT ont. ore: sy : 
«09 to of pie fssat' oat A nat 
fn20! edd nett veers! of . 
.biaaiwroo-ahev sil To mort aitwan 
‘sat yo beresam, Sh 7 Se 
“iebaved gytansyie ait 7 
stab adt .bentteb jesitllaliciods mM 

His px. | sshieeyinoe Fo, dsqnt ei rota : ners 
-trsqxe satan sit fom tee rol sobaanane Ye 
noftosysaon to taetofhiads Sat “ha gfdeT ak Hata oe 


orso byseaeD isan a oe i> 8 eel | 
meortequ aff mo Tealdysy tot Sonat ar 
ont 9907 gam bere eta a : wait be 
notsoiisenos afd Yo yetenamy fa (id 
vaneo ‘nots sorcsnnaie ie 

eweols oft to ces af one tosh 9 Ho 
wort barsowanaden wit? Yo oheey d\0 ett'bas teh 

smriad Fenotenemrb<aei apes oe ahtelige faotgor 


(1,3) (ond.9,m) + = 9 3 


4 4 . 


96 
tand, the slope of end-dump face,represents effect of closure material. 


The length of contraction can be expressed as 


(se) 
an aa F,(m,F tang) (5.2) 
Equation 5.2 can be written in an alternative form as 
we 
tamp, = SaeV Ply f(m,F tang) (5.3) 


A similar relationship for flow expansion can be assumed in the form 


L 
tanp, = TE BYE > F,(m,F tang) (5.4) 


Izbash (1970), based on experimental data, proposes a relationship for 


flow expansion in the form 
me 


which support the validity of equation 5.4. 

Based on the experimental data of this study, a plot of C. 
versus m is shown in Fig. 5.6. Because bo could not be measured pre- 
cisely, the scatter of the data is to be expected. The data differs 
from the hypothesis given by free streamline theory (Streeter, 1948) 
that C. assumes a minimum value of 0.611 for a two-dimensional slot. 

Wnekig. 5.7, PlOcs:.0F, tany, versus m and tany, versus 
[m“/10g(1/M) ] are shown. In Fig. 5.8, plots of tany, versus m and 


tany, versus [m“/1og(1/M)] are shown. Variation of F and tand are also 
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shown in these plots. 


As Figs. 5.6, 5.7 and 5.8 indicate, the correlations of tany, 


and tanpy with m is only fair, whereas the correlation of C. with m is 


good. The effect of variation of Froude number does not appear to be 


Significant, but the size of closure material appears to have a signi- 


ficant effect both on C. and tany,. Although, no definite conclusion 


can be drawn from these plots because of considerable scatter, the 


following observations can be made. 


(i) 


C. is primarily dependent upon m. The mean curve drawn through 
the data points in Fig. 5.6°can be used to determine Ce. Use 
of large and angular rocks can result in variation of C. 


upto + 10% from the mean value read from Fig. 5.6. 


(ii) Angle of contraction 1° as indicated in Fig. 5.7 varies 


between 5° and 25°. Use of large rocks would result in 
increase of V1: Precise prediction of its variation with 
other flow parameters is not attempted in this study, be- 


cause of limitations of the data. 


(iii) Angle of expansion bo as indicated in Fig. 5.8 varies be- 


tween 15° and 35°. Closure material size has a lesser in- 
fluence on Vo than Vy}: This is reasonable, as the expand- 
ing stream does not remember the effect of the contraction 
pattern. Froude number of approach flow seems to have a 
greater influence on bo than yy > It appears that an approach 
flow of higher Froude number, having a higher velocity of 


the jet at the vena-contracta, tends to draw the ambient fluid 
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in the eddy zone more quickly and to cause the jet to diffuse more 
rapidly. The process of diffusion and velocity distribu- 
tion downstream of the vena-contracta (Fig. 5.11) are similar 


to observations of Albertson et al (1950). 


The mean curves drawn through data points in Figs. 5.7b and 5.8b can 
only be considered as useful guides for determination of the flow 


profile at any stage of closure. 


5.2.3 Non-dimensional Plots of Contraction and Expansion Profiles 


From the live-stream configuration measured in each run, 
non-dimensional plots of flow contraction and expansion are shown in 
Fig. 5.9. (Fig. 5.5 is to be referred to as the key diagram). A non- 
dimensional plot of expansion due to Izbash and Lebedev (1961) is 
Shown on Fig. 5.9b for comparison. 

Figure 5.9 indicates that the profiles are similar, irre- 


spective of the contraction ratio and Froude number of flow. 


9.3 Velocity Distribution Through Constriction 
Figure 5.10 shows typical velocity distributions through the 


gap at the dam axis at increasing contraction ratios. Figure 5.11 
shows velocity distribution along the flow at one contraction ratio. 
For the runs shown in Figs. 5.10and 5.11 the uncontracted channel 
had a mean velocity of 1.15 fps, whereas the flow negotiating the 
closure dam attained a maximum point velocity of 1.81 to 2.05 fps. 
The interesting feature, however, is that although the contraction 


ratio increased from 0.30 to 0.75, the overall mean horizontal velocity 
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(b) FLOW EXPANSION 


FIG. 5.9 NON-DIMENSIONAL PROFILES OF THE LIVE STREAM BOUNDARY 
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in the contraction, as can be seen from Fig. 5.12, did not materially 
change, but remained around 1.50 fps, because of.considerable bed 
scour as shown in Figs. 5.10c, d and e. 

A pertinent observation that needs emphasis is that the 
velocity against the boundary as well as the stability of the material 
was relatively unaffected by scour at the dump line or downstream. 
Elaborating, let it be assumed that one size of material was found to 
attain the threshold condition at a particular contraction ratio. At 
that instant, the dumping was stopped. The scour progressed, thereby 
redistributing the velocity through the gap (decrease in mean velocity). 
After scour attained equilibrium, attempt to advance the dam further 
into the flow with the same material did not succeed significantly. 

The velocity distribution along typical sections such as XX 
and YY in Fig.5.10c were observed to be logarithmic. (The mean velocities 
for the contraction ratio of 0.65 shown were 1.83 and 1.57 fps re- 
spectively). The conclusive evidence was that the overall mean 


velocity through the contraction does not govern the stability of 


closure material at section XX. 


2-4 End-of-Closure Conditions 

At a contraction ratio generally in the range of 0.85 to 
0.90, the two wings of the closure dam join at the centre, after which 
the bed is no longer available for scour. Material instability then 
increases abnormally, because the afflux and flow velocity through the 


constriction suddenly increases abnormally (Fig. 5S.12)e 
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Izbash (1970) has observed this in actual closure jobs in 
the U.S.S.R. and designates this as the most critical stage of 
closure. According to him the flow attains the maximum specific 
power, which is the product of unit discharge through the gap and 
the head loss at the contraction. He recommends employing special 
blocks like tetrapods, or very heavy reinforced concrete blocks, or 
cribbed and garlanded boulders, to successfully close the flow. The 
design curves developed in this research program apply up to a con- 
traction ratio of 0.90 and cannot be used thereafter. It must be 
stressed here, however, that in actual practice at 90% closure,a 
very small fraction of flow really goes through the contraction, 
because most of it is averted through the alternate waterway. In 
rapid closures, discharge through the constriction will be less by 
the discharge going into storage. These aspects, combined with the 
small gap remaining, make the final closure a practical but difficult 


proposition. 


5.5 Stability of Closure Material 
5.5.1 Analysis of Stability Data by Critical Velocity Criteria 


In the preliminary runs, data at the threshold condition for 
each closure material were carefully collected to see whether instability 
could possibly be correlated with the experimentally observed maximum 
velocity occurring close to the boundary. A summary of three runs 
where threshold conditions were just attained is presented in Table 
5.1, which incorporates a comparison of the experimentally observed 


stable size and the size computed from the velocity formulas previously 
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TABLE 5.1] 
COMPARISON OF EXPERIMENTAL STABLE SIZE WITH 
VALUES GIVEN BY VELOCITY FORMULAS 


M Vmax D, D 
Contraction Observed Stable material size Experimentally 
Ratio Maximum in feet computed cor- Observed 

Velocity responding to velocity Stable 
Against Lees by Size in 

Boundary feet 
in feet 
per second Izbash Foc COrps Or 
(1970) Engineers 
(1970) 

] 2 3 4 5 
Ue25 1.80 0.02 0.025 0.021 
0.47 1.85 0.022 0203 0.059 
0.69 aS 0.023 0.038 0.082 


Hote: . For.al| the runs, 0.= 1.56 cfs..B = 7.50 ft... V.= lo 15 sips 
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109 
discussed herein. Correlation with mean velocity through the constric- 
tion (observed) has not been attempted, because this will indicate 
considerably smaller stable sizes. 

Good agreement at the low contraction ratio of 0.25 and 
poor agreement at higher contraction ratios of 0.47 and 0.69 between 
the experimental ly-observed stable size and the size computed accord- 
ing to velocity criteria can be noted in Table 5.1. This can be ex- 
plained as follows. 

At low contraction ratios, the increase in velocity and 
Shear stress on the boundary over the normal values is minimal. The 
flow behavior is still similar to that of a uniform flow in a pris- 
matic channel; whereas a higher contraction ratio results in increased 
flow curvature, flow velocity and acceleration close to the boundary 
and consequently results in a greater increase of shear stress (Ippen 
and Drinker, 1962), which demands a drastic increase in stable size. 
Experimental evidence (Table 5.1 ) substantiates the inadequacy of 
either the "mean velocity-through the constriction" or "velocity 
against the boundary" to govern the stability. 

In Appendix C, shear stress maxima occurring against the 
dumped material have been computed from measured velocity distribu- 
tion normal to the advancing boundary. Stable size of material has 
been computed from shear stress criteria, considering drag and lift 
(Christensen, 1971). Good agreement between experimentally observed 


stable size and computed stable size has been noted. 
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9.9.2 Proposed Design Curves 
(a) Effect of Turbulent Seepage 


It was observed that at high values of contraction ratio, 
seepage through the body of the closure dam became noteworthy because 
of increased head-differential across the dam. It was therefore con- 
sidered necessary to analyse the stability data after accounting for 
the seepage through the body of the dam. 

Detailed velocity measurements were taken at close intervals 
along the dump line from which, Qm, the discharge through the opening 
was computed. Assuming the discharge going into storage as negligible, 
Qo» the seepage discharge was computed as the difference of Q, the 
uncontracted flow discharge and Qn 

Because seepage reduces the discharge through the opening, 
its primary effect is to enhance the stability of the material on 
the end-dump face. Seepage, however, has a destabilising effect on 
the downstream slope (Terzagni and Peck 1948, Martin, 1970). 

In this study, the effect of seepage is accounted for by 
assuming the approach discharge to be Pe and not Q. It is, however, 
considered reasonable to disregard the effect of seepage on the 
stability of closure material, when Q. is less than 10% of Q. This 
is, of course, an arbitrary decision. 

It has been shown in section 3.1.5 that the essential non- 
dimensional parameters describing stability are F, D/h, mand n. It 
is obvious that a correct analysis of stability has to be based upon 


F and h values that correspond to Qn as the normal or approach discharge. 
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Detailed velocity measurements were taken in Series B, C, D, 

F and G. For runs in which seepage discharge exceeded 10%, the hydraulic 
particulars for both Q and Qn aS approach flow, are presented in Table 
2.2. B, the channel width is assumed the same for both cases. F and h 
corresponding to oe were obtained by referring to the stage-discharge 
curve of the channel for each series. 

The seepage discharge exceeded 10% of the normal, for closure 
material size 0.059 ft. and 0.082 ft. only. From a comparison of 
column 9 with column 12 and column 10 with 13 in Table 5.2 it is obvious 
that even a reduction of approach channel discharge (up to 30%) due to 
seepage does not result in appreciable (less than 10%) reduction of 
either V or hor F. A natural river also exhibits this trend, because 
reduction in discharge results in reduction of both stage and velocity; 
Froude number of flow is thereby reduced to a lesser extent. Because 
F and h are specifically considered in the analysis it is reasonable 
to expect that correlation of non-dimensional parameters based on 
normal flow data would yield realistic results except in the case of 
runs 34, 35 and 36. 

It needs to be mentioned that for all the runs in Table 
5.2, the contraction ratio was in excess of 0.70 and the efficiency 
of closure was as low as 0.70 or even lower. As will be discussed 
in section 5.5.2b, the design charts are developed considering efficiency 
up to 0.80. Although the data of runs in Table 5.2 are included in 
analysing the stability of closure material, the error introduced is 


considered minimal. 
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(b) Correlation of Essential Non-dimensional Parameters 

The summary of the experimental results incorporating the 
Salient data are presented in Table D.1. The maximum velocity Mage 
measured close to the boundary illustrates that its increase over 
increasing contraction ratios, is not as drastic as the increase in 
Stable size..demanded by the flow. 

The dimensionless variables n, m and D/h yielded by each 
run of series B, C and D are plotted in Fig. 5.13 for each of the 
three Froude numbers of 0.50, 0.29 and 0.10. The experimental evidence 
vindicates the validity of the four essential parameters provided by 
dimensional analysis. 

As the contraction ratio corresponding to the threshold con- 
dition could not be exactly determined experimentally, (and would also 
be a conservative estimate of material stability) a contraction ratio 
corresponding to an efficiency of 0.97 was arbitrarily considered as 
representating the onset of instability. However, the contraction ratio 
corresponding to an efficiency of 0.90 was accepted as the stable design 
condition forthe closure material, because loss of 10% of closure 
material will be economically commensurate with the gain in contraction 
ratio. Pariset and Hausser (1959) also accepted a loss of 10% of the 
closure material as marking the onset of an unstable condition. 

As has been mentioned earlier it is possible to obtain an 
improved closure of the waterway by sacrificing efficiency. This might 
well be the requirement in practical cases where only a partial closure 


is to be achieved by the economically available size. It can be seen 
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from Fig. 5.13b that for a D/h value of 0.204, an increase in contrac- 
tion ratio from 0.60 to 0.76 can be attained by allowing the efficiency 
to drop from 0595 to 0.60: 

As stated earlier, in section 3.1.5 "Dimensional Analysis", 
the omission of the term d/h is considered acceptable for practical 
design. 

Interpolating by numerical techniques, Fig. 5.13 has been 
transformed into Figs. 5.14, 5.15, and 5.16 as a three parameter plot 
of F, m and D/h corresponding to efficiencies of 0.97, 0.90 and 0.80. 
(An efficiency of 0.80 is acceptable in situations where availability 
of large size material is limited. Lower values may not be economical 
in practice.) Figure 5.15 is proposed as the design chart for stability 
of the closure material. Use of Figs. 5.14 and 5.16 can also be made 
as required. 

The data generally available for a closure job are Q, the 
river discharge, B, the river width, h, the normal depth and D, the 
closure material size. The information sought is b, the critical gap 
width. From the available information, F, the Froude number of approach 
flow and D/h can be computed and from Fig. 5.15, m, the stable con- 
traction ratio corresponding to an efficiency of 0.90, can be directly 
read out to yield b, the stable gap width. Alternatively, if the gap 
width b is specified, the contraction ratio m is known. Stable size, 


D, can then be evaluated by reading out D/h for known values of F and 


Mmeaga in trom Fig. 5.1 D. 
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FIG. 5.15 RELATIONSHIP BETWEEN D/h AND m FOR 
rn = 0.90 WITH F AS THE THIRD PARAMETER 
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FIG. 5.16 RELATIONSHIP BETWEEN D/h AND m FOR 
yn = 0.80 WITH F AS THE THIRD PARAMETER 
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(c) Effects of Neglected Parameters 
Effects of B/h Variation - The suggested design curves have 


been developed at one value of B (full flume width of 7.50 ft.) and 
one value of d (median size 1.20 mm). The effect of variation of B/h 
upon stability, as indicated by dimensional analysis, was studied by 
systematically varying B from 7.50 ft. to 4.6 ft. in series F and then 
to 2./5 ft. in series G. The experimental procedure for these runs 
was identical to those previously described. Only one Froude number 
of 0.29 was studied, utilising all the three sizes of closure material. 
Table 5.3 compares the observed contraction ratio with values read 
from Fig. 5.13. Although B/h varied in the range 28.8 to 10.6, no 
Significant difference in stability was observed over this range. 

Effects of d/h Variation - The effect of variation of d/h 
was specifically studied at a Froude number of 0.29 by varying the bed 
material size first to 0.60 mm in series J and then to 0.25 mm in series 
M. Other series in this study were run at Froude numbers of 0.50 and 
0.10 with 0.60 mm and 0.25 mm sand for collecting data specifically 
on scour. Because the stability data collected in those series were 
not very precise, effects of variation of d/h at F = 0.50 and 0.10, 
are not examined. 

To test the validity of equation 3.15, a semi-logarithmic 
plot of (p/h) (h/d) '/9 versus stable m (corresponding to efficiency of 
0.90 from data of series D, J and M) with F as a third parameter is 


shown in Fig. 5.17a. A linear correlation is indicated. 
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EFFECT LOR n VARIATION 


ye Experimental Values Contraction 
Run Be D Efficiency Contraction ratio m 
No. h h n ratio m From Figs. o. ls 
corresponding 
to efficiency 
n in col. 4 
] 2 3 4 5 6 
oy 17.70 0.081 0.95 0.435 0.45 
3) 17.470 0.081 0.65 0.65 0.64 
54 17.70 0.081 Orou 0.68 0.69 
95 1770 07227 0.98 0.54 0.50 
56 T7270 Oee27 0.80 0.69 Ol. #5 
57 17.20 0.027 0.60 0.74 0.78 
98 17.720 G3 0.95 0.70 0.72 
59 17.70 Oro) 0.65 0.75 0.82 
60 17.70 0.31 0.30 0.82 0.83 
6] 10.60 0.081 0.90 0.49 O51 
62 10.60 0.08] 0.90 Oey, 0.63 
63 10.60 0.227 0.98 0.44 0.50 
64 10.60 0.227 0.90 0.60 0.65 
65 10.60 0.227 0.50 0.68 Ogi, 
67 10.60 OS 0.90 0.69 0.74 
68 10.60 OE31 0.60 0.74 0.8] 


* - B/h = 28.80 for previous runs 
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FIG. 5.17 VERIFICATION OF EQUATIONS 3.15 AND 3.19 
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To dissociate the effects of h/d variation, a similar plot 
of D/h, m and F is shown in Fig. 5.17b. Although h/d varied in the 
range 47 to 320, for the Froude number of 0.29, the stability appears 
to be relatively unaffected. Table 5.4 compares the observed con- | 
traction ratio with values read from Fig. 5.13. The agreement is 
indeed good. 

It appears that for closure of alluvial channels the vari- 
ation of bed material size (consequently, the variation of bed scour) 
does not influence the stability of closure material significantly. 
The use of D/h in place of a term like (D/h) (ha) '/9 (equation 3.15) 
is justified in practical interest. Purely empirical correlations 
combining F and m to one variable are attempted in Figs. 5.16c and d. 
Because of the complex variation of the constants in equation 3.10, 
these linear correlations cannot be considered exact, but merely as 
guides. 

Figures 5.14 through 5.17 indicate that as the contraction 
ratio approaches unity, D/h values converge and the effect of the 
variation of Froude number is insignificant. 

It is concluded that the effects of B/h and d/h are minimal 
over the range tested. A possible explanation is that the contraction 
ratio and the geometry of the end-dump constriction can be considered 
primarily responsible for the distribution of shear stress (Ippen and 
Drinker, 1962, Replogle and Chow , 1966) and hence for the stability 
of the closure material. Variation of B/h has no significant effect 


on the contraction ratio. Variation of d/h results in variation of 
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TABLE 5.4 


EFFECT OF = VARIATION 


nn 


Experimental Values 


Run as D Contraction 
No Ai h Efficiency Contraction ration m 
: n ratio im TrOMe FIGs oO. 1S 

corresponding 
to efficiency 
Meine COL 

] 2 3 4 5 6 

US 132 0.081 0.98 0.39 0.38 

76 132 0.081 0.90 0.6] 0.68 

ld 132 0.31] 0.90 0.74 0.76 

79 132 0 Ziad. 0.90 0.65 0.66 

80 132 Ors) 0.80 0.80 0.81 

90 320 0.081 0.98 02.35 0288 

9] 320 0.081 0.90 0.50 Oe52 

92 320 0. 227, 0.90 0.65 0.64 

93 320 On) 0.90 Onse Ost 


*h/d = 45 and 66 for previous runs 
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the extent of bed scour, but as Fig. 5.1 indicates, the contraction 
ratio is not altered significantly due to bed scour. Hence the stability 
is relatively unaffected by variations of either B/h or d/h. 

(d) Effects of Rate of Dumping 

The rate of dumping though left out of the analysis, was varied 
widely in a few check runs to observe any effect on the stability. 
Until the onset of instability, the rate of dumping was observed to 
have no effect. After a particular closure material crossed the thres- 
hold condition, a higher dumping rate did not indicate significant 
improvement in efficienty, but merely resulted in the crest elongating 
along and not extending into the flow. In this study no quantitative 
assessment of the effect of either the rates of dumping or the variation 
of crest width has been made. It needs to be emphasized that a sub- 
stantially wider crest would make the head loss occur over a longer 
length and the closure material would gain some measure of stability. 
This aspect was taken advantage of by maintaining a 250 ft. wide crest 
for the Dalles dam (1965) in the final stages of closure. 

(e) Closure Tests from One Bank Only 

In series H, six verification runs were made at the three 
Froude numbers of 0.525, 0.29 and 0.1 to test the applicability of the 
design curves to a closure from one bank only. The contraction ratio- 
efficiency relationships were found to agree satisfactorily with values 
for dual closures, as can be seen from the comparison presented in 
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TABLE %5 5 


PARTICULARS OF CLOSURE TESTS FROM ONE BANK ONLY 


ee eee 


Experimental Values Contraction 
Run P D Efficiency Contraction Ratio m 
No. h n ratio, m from FIG sabals 
Corresponding 


to: thei nvGo war4 


1 2 3 a 5 6 
69 Oro *0-081 0.90 0.455 0.49 
70 On2e9e 202227 0.90 0.60 0.63 
7] Groce. blz 0-90 0.23 Oi 22 
ifs 0.525 0-327 0.90 0.40 0.46 
73 0.525 0.445 0.90 0.65 0.67 
74 QO310s 505056 0.90 0.74 0-73 


5.5.3 Additional Tests to Check Effects of Channel Cross-Sectional 
Shape 


It is well known that natural rivers rarely correspond to 
the rectangular or trapezoidal channel sections studied in the laboratory. 
It was therefore considered important to test closures in the laboratory 
on an undistorted model of a natural river. For this the Kitimat 
River, British Columbia, Canada, was modelled to 1:100 scale covering 
a stretch of over 3,000 ft. Closure was attempted with two sizes of 
material representing nominal diameters of 2.10 and 4.00 ft. at three 
discharges of 15,000, 30,000 and 50,000 cfs. Data for laboratory 


closure on a model of the outflow channel of Lake Athabasca, Alberta, 
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Canada (1971) are also available. A flow of 35,000 cfs was partially 
closed with stones representing nominal spherical diameter of 2.00 ft. 
The cross-sections of Kitimat River and the outflow channel 
of Lake Athabasca are shown in Fig. 5.18. As these rivers have a 
markedly deep central channel with side berms, the normal flow depth 
was substituted by the hydraulic mean depth for computing D/h values. 
A summary of the experimental data for these two cases are 
presented in Table 5.6. As can be seen, the agreement between the 
design value and the experimentally observed value of the critical gap 
width is quite satisfactory. It appears that the design curves can 
be applied to natural channels of irregular shape if hydraulic mean 


depth (R) is substituted for h, the normal flow depth. 


9.5.4 Checking of Design Curves 
(a) Against Actual Closures 


Comprehensive data of actual river closures or even partial 
closures by cofferdams are few due to lack of proper documentation. 
Many rivers in the U.S.S.R. (Linford, 1967) have been closed by either 
transverse-dump or end-dump methods, but exact data for end-dump dams 
pertaining to the threshold contraction ratio and the stable size are 
rarely available. Sandover (1971) quotes two cases, the first for a 
Russian river reported by Izbash, and the second for Mangla Cofferdam 
in Pakistan. A third case, Dalles dam on the Columbia river (1965), is 
also available, although it must be emphasized that the quoted value of 


critical gap width is rather approximate. A comparison of the critical 
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FIG. 5.18 IRREGULAR CROSS-SECTIONS TESTED 
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gap width from the design curve presented herein (Fig. 5.15) with that 
actually attained at the site is made in Table 5.7. 
(b) Against Data of Other Researchers 

Quantitative data on the stability of rockfill in end-dump 
closure of rigid-bed channels is available due to Sandover (1971). 
This is probably the only reported study on end-dump closures. Sand- 
over conducted the tests in a 10 ft. wide and 4 ft. deep flume, using 
two types of materials, rounded and angular (pebbles or chippings) in 
twelve different nominal diameters ranging from 0.3 inches for the 
pebbles and 0.3 to 8 inches for the chippings. Sandover's flume data 
are also utilized in Table 5.7. 

The agreement between the suggested design values and the 
actual river closure data as well as Sandover's data is indeed good. 
The actual river closure data are also plotted on Fig. 5.17d for 
illustration. 

However, all the closures verified were carried out on fixed- 
bed channels. The maximum contraction ratio attained in these cases 
was 0.89 for Mangla Cofferdam, but the Froude number was as low as 
0.035; under which condition for a mobile bed channel the bed scour 
would not be substantial, according to the indications of this study. 
The closure material on the firm bed channel would have then behaved 
similarly to that on the mobile bed channel, used herein to develop 


the design curves. 
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56 Concluding Remarks 


In this chapter the characteristics of flow through an end- 
dump constriction built into an alluvial channel and the stability of 
closure material were critically analysed based on the experimental 
data. Application of critical velocity criterion to determine the 
stable size was found to yield unrealistic results. 

A theoretical analysis of stability from shear stress 
criteria supported by dimensional analysis indicated the most important 
non-dimensional parameters to be the Froude number of approach flow, 
the contraction ratio, the ratio of nominal closure material size 
to depth of flow and the efficiency of closure. Experimental data 
corroborated the adequacy of these four parameters. Bed scour at the 
dump area, reflecting the effect of the bed material, was found to 
have a minor influence on the stability of closure material. Appendix 
C indicates satisfactory agreement between the observed stable size 
and the size computed from measured shear stress on the boundary. 

Curves correlating the essential parameters have been pre- 
sented as design curves for direct determination of stable size of 
closure material or of critical gap width at any stage of closure. 
These curves are not recommended for use at contraction ratios exceed- 
ing 0.90. 

A few actual closures and closures in the laboratory were 
analysed in the light of the proposed design curves. Good agreement 
between the actual critical gap widths and the predicted values is 


noted. 
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CHAPTER VI 
ANALYSIS OF DATA ON BED SCOUR AND BACKWATER RISE 


6.1 Introduction 

In this chapter, analyses of the experimental data on scour 
and backwater are made in the context of the theoretical developments 
made in Chapter III and the knowledge gained regarding the characteristics 
of flow through the constriction (section 5.2). 

Although it is accepted that characteristics of the bed 
material play a decisive role in the pattern and extent of scour and 
backwater rise, definite conclusion or design recommendations regard- 
jing local scour have been scanty, as discussed in Chapter II. The 
analysis of data in this study is made giving specific consideration 
to the size and the specific gravity of the bed material. 

Design curves correlating the essential non-dimensional 
parameters are presented. Conformity of the data of investigations 
on local scour by other researchers, with the proposed design curves, 
is examined. In such examinations, due consideration is given to the 
variation in the geometry and characteristics of the constriction from 


One situation to the other. 
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6.2 Maximum Scour Depth 
6e2 ie OCAtION= Ob Deepest Scour 


It has been shown earlier that at the vena-contracta the 
contracting flow attains the maximum velocity close to the live stream 
boundary and the flow velocity decreases towards the centre of the 
channel (Fig. 5.11). Two isolated elliptical scour holes with some 
general scour of the bed in the centre generally occur at low con- 
traction ratios. As the contraction ratio increases the scour holes 
move to the centre and become deeper. Figures 6.1 and 6.2 illustrate 
the typical scour pattern at different degrees of contraction at a 
Froude number of 0.50, for the bed material of 1.2 mm. The same trend 
was observed at the Froude numbers of 0.29 and 0.10 and for the bed 
materials of 0.6 mm and 0.25 mm as well. From the definite trend in 
the location of the deepest scour close to the live-stream boundary, 
at the vena-contracta both for the clear water and the sediment- 
transporting flows, it is obvious that the increased bed shear stress 
due to the contracted flow is primarily responsible for maximum scour 
at the end-dump constrictions. 

Figures 5.7 through 5.9 which define the live-stream con- 


figuration, can be used to determine the location of deepest scour. 


6.2.2 Clear Water Case 
The data for the runs 7 through 85 and 94 through 96 corr- 
esponding to the clear water case are presented in Table D.1. 


The analysis of the data is based upon Eq. 3.31 and Eq. 3.48 
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A semilogarithmic plot of (h +d,)/h versus m with F as the third para- 
meter is shown in Fig. 6.3a. The plot indicates a definite correlation 
between (h + d.)/h, mand F. On this figure, the variation of h/d is 
Shown to illustrate that the sediment size has a definite influence 
on the maximum scour depth. 

(a) Effect of Sediment Size 


In Fig. 6.3b a plot of {(h + d.)/h}(d/h)°° 


| versus m shows 
the data to coalesce at each Froude number which illustrates the validity 
of Eq. 3.28. The index 0.29 indicates that at a particular Froude 
number of the flow, the variation of the sediment size in the practical 
range from medium to coarse sand (say from 0.2 mm to 1.5 mm, 750% 
increase) would result in a decrease of the non-dimensional scour depth 
by only 75%. However the bed forms which would be more pronounced for 
the fine sand can cause variation of scour depth over time. 
It has been recognized that the critical tractive stress 
for dune beds is higher than that for plane beds (ASCE Task Committee 
on Sedimentation, (1966)). As such Eq. 3.29, derived under the assump- 
tion of an invariant critical tractive stress related to size only, 
would not be completely realistic for flows over different bed forms. 
At low values of particle Reynolds number, Eq. 3.29 is not applicable. 
These considerations and factors such as the shape of sediment can 
cause the index of 0.29 to be as low as 0.16 as suggested by Neill (1962). 
Figure 6.4 shows the non-dimensional plot of experimental 
data, based on Eq. 3.31. The theoretical analysis has thus a satis- 


factory agreement with the experimental data. The combined parameter 
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FIG. 6.3 RELATIONSHIPS FOR MAXIMUM SCOUR DEPTH IN CLEAR-WATER FLOW 
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ADDITIONAL LEGEND FOR CHECK RUNS 


F inmm h/d 
(a)FOR PEBBLES @ 0.70 66 10 
@ .050=do 20 
Ay 0.50. dona l0 
© 050 do 15 
© 036 do 13 
m 029 do 20 
(b) FOR COAL e 015 16 57 : 
e@ 0.10 


FIG. 6.4 VERIFICATION OF EQUATION 3.31 FOR MAXIMUM 
DEPTH OF SCOUR IN CLEAR-WATER FLOW 
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{(h + d_)/h}(a/n)?*4? seems to represent the effect of sediment size 
more adequately than Cas the average drag coefficient of the sediment 
as proposed by Garde et al (1961). It was, however, considered de- 
Sirable to check the validity of Eq. 3.31 with coarse bed material 
and light weight bed material for general applicability. 
(b) Check of Data for Large Sediments 
It was considered desirable to check Eq. 3.31 with pebbles 
of 6.6 mm (closure material B) as the bed material. Table 6.1 sum- 
marizes the data for the check runs. The data is also plotted on 
Fig. 6.4 for illustration. It is seen that the agreement is satis- 
factory for h/d values as low as 10. 
(c) Check of Data for Light-Weight Bed Material 
Crushed coal (anthracite) of median size 1.6 mm and specific 
gravity of 1.26 was used as bed material in a few check runs. The 
Summary of the experimental data is also incorporated in Table 6.1. 
The data is also plotted on Fig. 6.4. Satisfactory agreement with 


Eq. 3.0!1- 1S nosed, 
(d) Maximum Scour Depth in Alluvial Channels 


As the experimental data agrees with equation 3.28 which in- 
corporates the essential dimensionless parameter given by Eq. 3.48, 
Fig. 6.4 can be used to determine the maximum depth of scour. To 
obviate tedius repetitive mathematical computations for determination 
of (d. + h)/h for alluvial channels from Fig. 6.4, convenient charts 
covering Froude numbers in the range 0.05 to 0.60 are presented in 


PiGs 8020 =O DeG sediments of specific gravity 2.65. 
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(e) Evaluation of Indices J, and Jy in ce /t-m" Relationship 
max 


Based on Eq. 3.31, the slope of the straight line fitting 
the entire data in Fig. 6.4 is 0.43 j3 and the intercept on the ordinate 
is 0.43 Jy: Thus J3 and Jq are evaluated as 2.1 and Sais respectively. 
Awazu (1967) measured tote at the centre line of the flume and evaluated 
J3 as 1.40 and J4 as - 0.021. Because velocity and shear at A is greater 
than that at B (centre of the gap) in Fig. 2.2 , the results of this 


Study are compatible with Awazu: however, it is seen that Ne 
max 


Eq. 3.26 is primarily dependent upon the contraction ratio and is largely 
independent of the Froude number of the flow. This was also observed 
in drag measurements on end-dump face (Appendix C). The other interesting 


observation is that up to a contraction ratio 0.10, [3 , does not 
max 


exceed unity. The exponential form of Eq. 3.26 can be assumed reason- 


ably valid in the contraction ratio range of 0.15 to 0.85. 


6.2.3 Sediment-Transporting Flow 
In the runs 86 through 93 and 97 through 100 the flow was 


transporting sediment. The experimental data are presented in Table 
Daly, 

Based on Eq. 3.43, a semilogarithmic plot of (h + d.)/h versus 
m with F as third parameter is shown in Fig. 6.6a. The data shows a 
systematic variation with F. Although the effect of d/h has not been 
studied specifically, it appears to be insignificant in the range 
tested. In fact for F remaining nearly constant (0.50 to 0.45), with 


the increase in h/d from 91 to 232, (h + d.)/h decreases, instead of 
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FIG. 6.6 DESIGN CURVE FOR MAXIMUM SCOUR DEPTH 
IN SEDIMENT-TRANSPORTING FLOW 
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144 
increasing as would have been the case with clear water flow. It 
needs to be mentioned that the experimental procedure for these runs 
was deficient, because of the non-circulating nature of the flow. 
Manual feeding of the sediment at the upstream as was actually done, 
can not be considered to maintain absolutely uniform flow and sedi- 
ment discharge from.upstream. In view of these a truly rigorous 
analysis is not attempted. 


In Fig. 6.6b, a combined parameter (h + d.)/h ey 


0.85 Occursfinseqmic. oe 


plotted against m. (For clear water flow F 
The correlation being satisfactory, this figure is proposed as the 


design curve for sediment-transporting flow. 


6.3 Time History of Scour 
O-o. I Characteristic time sof Scour 

It has been recognized that the rate of scour is influenced 
by the sediment properties. Rouse (1939), Laursen (1953), Doddiah, 
Albertson and Thomas (1953), Feld (1968), Govindrao and Sharma (1967) 's 
analyses of scour due to jets indicated V/w as a significant para- 
meter describing the rate of scour. Carsten (1966) indicated that Vt/h 
in combination with d/h describe the rate adequately where h is a 
characteristic length of the system. Many investigators (Ahmed, 1953, 
Breusers, 1967, Watkins, 1969) have verified that the depth of scour, 
after the initial transitional period (when the rate is extremely 
rapid) has a linear semi logarithmic relationship with time. However 


a definite relationship of the scour rate and the flow, fluid, characteristics 


and the geometry of the obstruction has still not been obtained, because 
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of the complexity of the problem. Dimensional analysis and experi- 
mental data have been used with judgement and as Breusers (1967) 
indicated, exact scaling of the model data to the prototype still 
remains unsolved particularly for fine sediments. Breusser indicated 
that the time scale of scour can be expressed as a function of the 
length scale, sediment size scale (d prototype/d model), density scale 
and velocity scale [(V-V,.) prototype/ (V-V,.)mode1]. The functional 
relationship will vary from one situation to the other. However a 
similarity in the scour hole configuration does occur (Laursen, 1952, 


Breusers, 1967) and can be expressed as 


Se (ee) (Cap 


in which d, and x represent the characteristic depth and length of 
the scour hole at time t. The determination of the equilibrium time 
te for any scouring process is the most difficult problem. 

Bed profiles recorded in this study are analysed to develop 
relationships for the rate of scour for the clear water case only. A 
few typical plots of depth of scour versus time are shown in Fig. 6.7. 
Semilogarithmic plotting of the same data against the non-dimensional 
parameter Vt/h (Eq. 3.47)showasystematic trend (Fig. 6.8). This 
indicates that with F, m, d/h and SI TekG kept constant, (h + d,)/h 
should be uniquely related to Vt/h. Other researchers have made 
similar assumptions and have obtained compatible results. For example, 


Laursen (1963) postulated that d, represents the hypothetical maximum 
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scour depth for a coarser sediment to develop an expression for the 
rate of scour. Govind Rao and Sharma (1967) also accepted a similar 


hypothesis and proposed for scour due to a vertical jet 


0.4 
tan7! 9.4 (Wty LT (6.2) 


n 
in which T and T' are constants. A similar expression for scour due to 
a horizontal jet was also developed empirically. T andT' are evaluated 
from the boundary conditions of d,. = 0 at t = 0 and d, = d. at t = ~, 
Equation 6.2 satisfies the scour relationship over the entire period 
and has the advantage of yielding the maximum scour depth directly. 
Based on the foregoing considerations and Eq. 3.47, a non- 
dimensional form of the transient-scour depth can be deduced from Eq. 
3.31, with values of 53 = 2.1 and J4 = - 0.17, as 
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Particulars of representative runs, covering three bed materials at 
three Froude numbers, chosen for verification of Eq. 6.3 are presented 


in Table 6.2. 
Bplot of Eq..6/3 in Fig. 6.9 Indicates asdefinitertrend 


though showing a wide scatter. This is understandable because the 
scouring process is certainly not as simplified as assumed in this 


analysis. Figure 6.9 is not proposed as a design curve because of its 
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TABLE 6.2 
PARTICULARS OF RUNS ANALYSED FOR RATE OF SCOUR 


Run Q B V h d F b 
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FIG. 6.9 VERIFICATION OF EQUATION 6.3 FOR SCOUR 
DEPTHS DURING THE ACTIVE PHASE 
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limitations. It can only be used as a guide to determine d, approxi- 


mately at any time t. 


6.3.2 Dimilarity in Scour 


From a mean curve plotted through the data points in Fig. 
6295 tgo> the time to attain 90% of the maximum scoured depth was read 
out for the runs in Table 6.2. A dimensionless plot of t/t versus 
(d, + h)/0.90(d. + h) for all runs (Fig. 6.10) show that scour pattern 
for a particular geometry (end-dump constriction in this case) is 
Similar irrespective of the bed material size. This was also observed 
by Breusers (1967). An important conclusion from Fig. 6.10 is that 
before 0.1 tg is elapsed, 60% of the maximum scour occurs, which 
illustrates the necessity of execution of a closure job with minimum 
interruption. For practical purposes, an interruption of a few hours 
on the prototype, would be expected to result in development of full 
Scour in the dump area. 

A typical plot nondimensionalised with respect to the length 
and depth of the scour hole is shown in Fig. 6.lla. Similar plots 
were obtained for all runs. The different profiles of the scour 
hole with respect to time fall on one curve at all values of non- 
dimensional time as shown in Fig. 6.11b. This was also the conclusion 
reached by Laursen (1952). From the pattern of the scour profiles 
(Figs. 6.1 and 6.2), it would be reasonably correct to expect the 
deepest scour to occur at the vena-contracta section, close to the 


live-stream boundary. The contraction of the live stream and the 
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location of the vena-contracta can be assessed from Figs. 5.7 and 5.9. 
Thus the location of the deepest scour hole can be ascertained with 
reasonable accuracy at any stage of closure. The configuration of 

the upstream end of the scour hole can then be arrived at by assuming 
the upstream slope of the scour hole equal to the angle of repose of 
the bed-sediment (Fig. 6.1la). The side slopes were observed inclined 
at the angle of repose in this study. The downstream slope, which 
both Carsten (1966) and Laursen (1963) assume to be the same as angle 


of repose was however observed in:this study to be much flatter. 


6.4 Check of Design Curves Against Previous Studies 
6.4.1 General 

Investigation on scour at end-dump constrictions do not 
appear to have been reported. Relevant studies on local scour at 
constrictions are chosen to conduct checks. 

Two types of checks are made to verify tne agreement be- 
tween the data leading to the design curves developed in this study 
and the recommendations or data of previous studies. Pata of clear 
water runs in this study, being considered more reliable are analysed 
only. 

(1) The data of this study in respect of maximum scour depths 
are compared with 

(a) Laursen's (1960,1963) design curves for i) ioca! scour 

at an isolated abutment, and ii) general scour within 

o posing abutments. 

(b) Komura's (1966) design curves for general scour at a 


long contraction. 
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(2) The maximum scour depths measured at 
(a) vertical-edge spur-dikes (Garde et al, 1961) and at 
(b) bridge constrictions (Liu et al, 1961) are compared 


with this study. 


These studies did not distinguish local scour from general 
scour and design curves were proposed for the particular geometry 


of constriction studied. 


6.4.2 Check Against Laursen's and Komura's Analyses 
(a) Isolated Scour Hole (due to Laursen, 1963) 


A comparison of the observed d/h for a few typical runs of 
this study having m less than 0.65, (considered equivalent to an isolated 
abutment) is made in Fig. 6.12 with Laursen's curves (shown in Fig. 2.3). 
It can be seen that the maximum depths of scour obtained in this study 
are consistently lower (1/5 to 1/2) than the values predicted by 
Laursen. Thus, the maximum depth of an isolated scour hole at an end- 
dump dam can not be obtained from Laursen's analysis for isolated 
abutments. 

(b) Overlapping Scour Hole (due to Laursen, 1960 and 1963, and 
Komura, 1966) 

The maximum depths of scour pertaining to a few typical runs 
of this study, having m in excess of 0.65 are considered appropriate 
to be analysed as data of scour in a long contraction. These are 


plotted in Fig. 6.13a for check against Laursen's curve and are plotted 
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in Fig. 6.13b for check against Komura's curve. The contracted width 
b. (measured in this study) is considered equivalent to the contracted 
width in a long contraction. The agreement of the data both with 
Laursen's and Komura's curves is quite striking which indicates that 
an end-dump constriction behaves as a long contraction at higher con- 


traction ratios. 


6.4.3 Check of Garde et al's and Liu et al's Data 

The curve in Fig. 6.3b is redrawn on Fig. 6.14 combining the 
non-dimensional parameters (d.th)/h, F and d/h to one parameter based 
on equation 3.31. 

The observed maximum scour depths in a few typical runs from 
Garde et al's studies on bluff vertical spur-dikes and from Liu et al's 
(1961) studies on spill-through abutments are plotted on Fig. 6.14 for 
comparison with this study. The particulars of the runs chosen are 
detailed in Table 6.3 in which the observed and computed depths are 
compared. 

The observed scour depths are consistently higher than the 
values predicted by this study, although the difference is less at 
higher contraction ratios. A similar trend was observed in the com- 
parison with Laursen's curves for isolated abutments. Because an 
end-dump dam has a better stream-lined shape than a bluff, rigid 
abutment or spur-dike, a reduction in the maximum depth of scour is 
reasonable to expect. The results of this study, therefore, appear, 


compatible with the previous studies. 
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TABLE 6.3 


COMPARISON OF GARDE et al's (1961) DATA AND 
LIU et al's DATA WITH THIS STUDY 


Sn eee Use eee oS eee eI Dee ee De 


d F m 

Run Size of Froude Contraction (h+d_)/h (h+d_)/h 
No. bed Number Ratio : ae 

Material of normal (Observed) Determined 

in mm Flow from Fig. 

6.14 
] 2 3 4 5 6 

Garde, Subramanya and Nambudripad (1961, 1963) 
Vertical edge spur-dike 
D-9 0.29 0.14 E47 2.08 Vet2 
D-10 0.29 0.284 0.47 Seco 2.60 
D-12 0.29 0.23 0.47 easy 2820 
A-3 0.29 On 65 0.33 ei Talo 
A-10 0.29 0.12 0%33 1238 02 
E-] 0.45 0.148 0.165 1.29 No scour 
E-3 0.45 Ov225 0.165 be/9 No scour 
E-7 0.45 Gaciic 0.165 2.00 1.08 
F-2 1.00 0.172 0.165 Lead No scour 
F-3 1.00 0.24 0.165 1.68 No scour 
F-6 1.00 0.318 0.165 ea a | 1305 
G-1 ae vhs) 0.182 0.165 1202 No scour 
G-2 225 0.408 0.165 1.61 No scour 
G-3 ta iape 0.476 Oat65 1.80 1.06 


Liu, Chang and Skinner (1961) 
Spill-through model 


pe 


146 0.56 0.242 0.47 cm 193s 
148 0.56 Once 0.47 SENOS aS ane) 
180 0.56 0.166 Cor ots No scour 
182 0.56 Vacs) (ORE: 2.63 Lado 
203 0.56 0.16 0.50 Ta90 2.03 
206 G56 0.394 0.50 Sater eee) 
210 0.56 0.664 0.50 5.40 4.60 
232 0.56 0.203 Ose leis ERAS 
234 Jeo 0.378 0.37 3.00 2.80 
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6.5 Proposed Design Curves 

In view of the good agreement of the experimental data with 
the theoretical analysis both for the clear-water and sediment-trans- 
porting flows and satisfactory checks with previous studies, Figs. 
6.4, 6.5 and 6.6 are proposed as design curves for determining maxi- 
mum scour depth. 

The experimental data indicates F, m, P./ Ps d/h and 
(d, + h)/h to be the essential inter-related non-dimensional variables 


describing scour. 


6.6 Maximum Backwater Rise 
6.6.] Location of Maximum Backwater | 

The basic data is presented in Table D.1. Typical water 
surface profiles through the constriction have been shown in Fig. 5.2 
which indicates that although the section of maximum backwater is 
some distance upstream of the dump line, the level in the upstream 
deadwater zone at A (Fig. 5.2) very nearly corresponds to the maxi- 
mum backwater level. (H is measured accordingly in contrast to h, 
h Fig. 1. 2b) which are measured along the center line 


h, and h 


me 2 3 ( 
of flow. Difference between H and hy is not appreciable). 

Liu et al (1957) and Biery and Delleur (1962) have corre- 
lated the length to maximum backwater Lio (Fig. 1.2), gap width b, 
Froude number F and opening ratio M. It is generally seen that 


the non-dimensional term Ly2/b decreases with the increase in Froude 


number. It also decreases with the increase in the opening ratio. 
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For practical purposes L,,/b can be considered to vary in the range 
0.5 to 2.5. Carter and Tracy (1955) observed that in general Lyo/b 
can be taken as unity. 

From the observation in this study, it can be stated that 
exactly locating the maximum backwater section is not possible, 
because of the extremely low variation of slope in this region. 
From practical considerations, the water level just upstream of the 
dam should be considered equal to the maximum backwater level for 
fixing the crest level of the dam, with suitable freeboard. The 
discharge capacity of any diversion arrangement, which is generally 
sited a short distance upstream of the dam, should be arrived at 


considering the maximum backwater level as the headwater level. 


6.6.2 Variation of (H/h)> with F°(1/M* - 1) 

A plot of (H/h) versus m with F as the third variable is 
Shown in Fig. 6.15. The h/d values are noted on this figure. For 
clarity the data points for the different sizes of closure material 
are not distinguished but it is seen that at the Froude numbers of 
0.29 and 0.50 for the same h/d values the maximum backwater rise is 
slightly reduced with the larger closure material because of in- 
creased seepage through the body of the dam. 

Figure 6.16 shows a plot of (H/h)> versus Fo(1/Me - 1). 
The effect of h/d variation does not seem to be significant, both 
for the clearwater and sediment-transporting cases. In fact it is 


seen that an increase in h/d which results in an increase in the maximum 
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FIG. 6.15 VARIATION OF BACKWATER RATIO H/h, WITH CONTRACTION RATIO m 
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depth of scour results in a slight increase of the backwater level. A 
valid explanation regarding this has not been found. On the other hand 
the relative constancy of the maximum backwater level with varying h/d 
can be explained as follows. For a strictly free fall or flow down 
a steep slope the flow would pass through the critical before the out- 
fall and the outfall functions as a control. Although the flow depth 
does not become critical so long as the bed is available for scour, the 
flow approaches critical condition with increase in contraction ratio. 
It is therefore possible that the upstream depth is primarily related 
to F and m at higher contraction ratios and an increase in depth of 
scour downstream does not significantly lower the water level upstream. 
But the maximum backwater rise for a mobile bed is certainly lower 
than that for a rigid bed, as the curve due to Liu, Bradley and Plate 
for spill through abutment (comparable to end-dump constriction) 


shown on Fig. 6.16 indicates. 


6.6.5 8LVariatvion of C with F and m 


Assuming equation 2.25 to be valid, the discharge coefficient Ca 
for all runs were computed. It is apparent that a correlation of Ce 
F and M exists for flow through constriction on mobile beds. Omitting 
the data points, Fig. 6.17 is presented as a plot of C versus F with 
mas third parameter. Here again the effect of h/d variation has not 
been considered separately because the data does not indicate any 
systematic trend. On Fig. 6.17, the corresponding plot for flow on a 


rigid bed due to Valentine is shown for comparison. Figure 6.17 indi- 
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cates that as the flow passes the dump line and continues down the 


scour hole, the coefficient is considerably increased due to suction. 


6.6.4 Design Curves 


Because the effect of particle size on the maximum back- 
water rise is not conclusively studied, the curves in Figs. 6.16 and 
and 6.17 can not be considered as being precise but rather as useful 
guides in practical cases. However, the design curves are considered 
to have better applicability than the tentative recommendations of 
Liu, Chang and Skinner, who experimented with 0.56 mm sand only and 
stated "for the sediment involved in the tests, the water-surface drop 
across the embankment is approximately 60% of that in a rigid channel. 


It is expected, however, that this ratio varies with the bed material". 


6.7 Concluding Remarks 


In this chapter, the experimental data on scour and back- 
water have been analysed. 

Assuming an exponential increase of the ratio of the maximum 
bed shear stress in the constricted channel to the boundary shear 


stress of the normal flow with the contraction ratio cae as with 
ma X 


m), an analysis for the maximum depth of scour was made in Chapter III, 
(section 3.2 ). The experimental data on scour supports the theoretical 
analysis well. Design curves for auc the maximum scour depth 
for the clear water and sediment-transporting flows have been proposed. 
Similarity in scour profiles both in respect to space and time were 


observed. Tentative design curves for the time history of scour have 
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been proposed. The size of the bed material appear to have signifi- 
cant influence both on the depth and rate of scour for clearwater 
flows. For sediment-transporting flows bed material size has less 
influence on the depth of scour. The most Significant inter-related 
non-dimensional variables describing scour are the Froude number of the 
approach flow F, the contraction ratio m, the ratio of the bed ma- 
terial size to depth of flow d/h, the specific gravity ratio Petey 
and the ratio of scoured depth to original depth (d.t+h)/h). 

Maximum backwater rise due to constriction on an alluvial 
bed was observed to be less than that for a rigid bed, as obtained 
in the studies of Liu, Bradley and Plate (1957) and Valentine (1958). 
The depth of scour does not appear to significantly affect the back- 
water rise (probably beyond a certain depth of scour), as does occur 
in the case of flow over a free fall. Definite conclusion regarding 
the effect of bed material size on backwater rise could not be 
drawn. Froude number of the approach flow F, the contraction ratio 
M, the coefficient of discharge C. (assuming a weir flow through 
the constriction) and the non-dimensional backwater rise H/h were 


observed to be the significant inter-related variables. 
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CHAPTER VII 


SUMMARY , CONCLUSIONS AND RECOMMENDATIONS 


7.1 Summary and Conclusions 
7.1.1 Introductory Remarks 


The present study has been undertaken to develop comprehen- 
Sive hydraulic design procedures for end-dump closure of alluvial 
Channels at right angles to the flow. The primary problem in a river 
closure is the determination of the stable size of material (generally 
rockfill) at all stages of closure. The literature review indicated 
that the stable size is generally computed from the maximum velocity 
through the constriction, applying the "velocity-size" correlations 
developed for transverse-dump closures by Izbash (1936). This pro- 
cedure was considered to be of questionable accuracy, because the flow 
against the sloping face of an end-dump dam is curvilinear and rapidly 
accelerating in contrast to the rectilinear flow over the level bed 
of a transverse-dump dam. 

Specific studies on end-dump closure of rigid-bed channels 
leading to design formulae for the stable size have been reported by 
Pariset and Hausser (1959), and Sandover (1971). These formulae do 
not relate the stable size to velocity, but contain non-dimensional 
numbers describing the flow, fluid, and closure material. Applica- 


tion of Pariset and Hausser's formula requires predetermination of 
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the backwater rise and head loss at the constriction and the authors 
have not suggested any method of evaluating them. Sandover's pro- - 
cedure involves a trial and error solution. As such both procedures 
are impractical as a means of directly determining stable size during 
closure of alluvial channels. A systematic study on end-dump closure 
of alluvial channels does not appear to have been reported. 

The literature review also revealed the lack of systematic 
investigation or specific design recommendations in respect to the 
associated phenomena such as the flow pattern, bed scour and the back- 
water rise at end-dump constrictions on alluvial channels. 

An attempt has been made in the present study to formulate 
relationships for closure material stability, scour, and backwater, 
involving all of the possibly relevant variables. These relation- 
ships have been evaluated and design recommendations have been made 
on the basis of experimental data from laboratory tests consisting 
of over 120 runs. Three sizes of closure material (pebbles and 
chippings, 6.60 mn, 18 mm and 25 mm) and three sizes of bed material 
(sand, 1.20 mm, 0.60 mm and 0.25 mm) were utilized in the tests, which 
were carried out over three Froude numbers of the approach flow (0.50, 


0.29 and 0.10). 


(ol.ePeotabillaicy of Closure Material 


At the outset of the experimental programme, the flow pattern 


was observed and the velocity distribution was obtained through an end- 
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dump constriction at different contraction ratios. The experimentally 
observed stable size was found to be consistently and substantially 
larger than that computed by correlation with the maximum velocity 
observed against the closure material. It was clearly demonstrated 
that neither the mean velocity through the constriction nor the maxi- 
mum velocity against the closure material adequately describe the 
hydrodynamic effect of the curvilinear, accelerating flow. 

A numerical measure termed “efficiency of closure material" 
has been introduced to denote the volume of closure material remaining 
Stable within the limits of the dump line as a percentage of the 
quantity dumped into the flow at any stage of closure. 

In Chapter III, an analysis of stability was done by assuming 
the closure material to become entrained after the hydrodynamic forces 
on the end-dump face exceed a critical value. Christensen's (1971) 
expression which includes the effects of particle lift, and bank slope, 
was used to determine the critical shear stress. A dimensionsal 
analysis was made of all the pertinent variables characterizing the 
flow, fluid, closure material and bed material. These analyses in- 
dicated the essential non-dimensional parameters to be the Froude 
number of the approach flow F, the contraction ratio m, the ratio of 
the closure material size to the approach flow depth D/h, and the 


efficiency of closure n. The form of the function which was subjected 


to empirical evaluation was 


0 = f,(F, m, D/h, n) (7.1) 
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In Chapter V, data on the flow characteristics and material 
Stability have been analysed. Design curves (Figs. 5.5 through 5.9) 
to arrive at the live stream boundary covering the flow contraction 
and expansion have been presented. The contraction ratio was ob- 
served to be the most important parameter describing the flow configu- 
ration. Variation of bed scour at the dump line at one contraction 
ratio did not appear to have appreciable effect on the stability of 
the closure material. In other words the bed material size, though 
influencing the depth of scour considerably, did not influence the 
Stability on the end-dump face significantly. This was checked by 
varying h/d in the range 45-500. Similar variation of B/h in the 
range 10-29 did not influence the closure material stability signifi- 
cantly. 

Design recommendations in the form of curves (Figs. 5.14 
through 5.16) correlating F, m, D/h and n have been presented from 
which the stable size of closure material of standard shape and 
density can be determined directly. These curves are not recommended 
for use beyond a contraction ratio of 0.90. An efficiency of 0.90 
was taken to adequately represent the stable design condition. The 
applicability of the design curves to irregular cross-sections of 
natural rivers has been tested by assuming the mean hydraulic radius 
to be representative of the normal flow depth. A few actual proto- 
type closure data and the experimental data of Sandover have been 
compared with the recommendations of this study and the agreement 


was satisfactory. Satisfactory applicability of the design curves 
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to closure from one bank only was also noted. 

From measured velocity distributions normal to the end-dump 
face, the maximum shear stress Lotan on the boundary was computed from 
boundary layer theory. The shear maxima was observed to follow the 
maximum velocity flament. An exponential increase of "t ioe! 
with "m" was found. (Appendix C). The experimental Tee he was 
observed to agree reasonably with the size computed to be stable against 


Ty) considering both drag and lift. For instability on the end- 
max 


dump face Shields dimensionless shear stress Tt, was found to 


assume a value of 0.02. 


721.0 ~bed Scour 

The scour patterns observed in this study confirmed that the 
bed scour at an end-dump constriction can be considered to be a com- 
bination of local scour (occurring at obstructions) and general scour 
(occurring in channel contraction). The self-formed end-dump face, 
however, stream-lined the flow, so that the deepest scour was observed 
to always occur downstream of the dump line. A reasonably two-dimensional 
flow was observed at re vena-contracta, where the maximum flow velocity 
and the associated maximum boundary shear resulted in the scour hole 
attaining its maximum depth. The term “maximum scour depth" represented 
by d. has been used herein to denote the limiting depth of scour below 
the original bed in the case of clear water flow and the equilibrium 
depth of scour below the original bed in the case of sediment-transporting 


flow. 


The previous studies on scour at constrictions such as spur- 
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dikes, and bridge abutments and at long contractions have been reviewed 
in Chapter III. The information was found to be deficient and incon- 
clusive as to the effects of both the bed material size and density. 
Simplified analyses for arriving at the maximum depth of scour were 
made in this study both for the clear-water case and the sediment- 
transporting case by relating the rate of scour to the increased bed 
Shear stress at the vena-contracta under an assumption of two-dimensional 
boundary layer flow. Maximum scour depth in clear-water flow was 
assumed to have been attained when the boundary shear in the scour 
hole is reduced to the critical value for the bed sediment, whereas 
in sediment-transporting flow it is attained when the rate of trans- 
port into and out of the scour hole equalizes. Dimensional analysis 
indicated the form of relationship for clear water flow (equation 3.31) 


as 
(htd.)/h = f5(F.m,h/d,ow/Ao,) (7.2) 
and for sediment-transporting flow (equation 3.43) 
(h+d.)/h = f3(F,m) (793) 
Experimental investigations led to the correlation of the 
dimensionless terms in equations 7.2 and 7.3 and design curves (Figs 6.5 


6.6 ) are presented for determining the maximum depth of scour. The 


applicability of the design curves has been tested for light-weight 
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material (coal) and 6.60 mm size pebbles as bed material. Data of 
other researches on scour at spur-dikes, bridge abutments and long con- 
tractions have been compared with the results of this study and are 
compatible with the test results obtained here. 

The rate of scour during the active phase was observed to 
be logarithmic with time and tentative design curves have been pre- 
sented to determine the transient scour depths approximately. Because 
the maximum depth of scour is reached asymptotically, time to attain 
90% of the scoured depth, tgo> has been considered as a characteristic 
time and an important observation was that within an interval of 0.1 
tg from the commencement of scour more than 60% of the maximum scour 
occurred. A similarity in the scour hole configurations with respect 


to both space and time were observed. 


7.1.4 Backwater Rise 
Based on the analysis of backwater rise by previous investi- 
gators for constrictions in rigid-bed channels, two different forms 


of backwater rise equations were framed for experimental investigation. 


(H/h)° fF - 1)) (7.4) 


and 


Wed Ue pte (7.5) 


in which Ci. the coefficient of discharge for the flow through con- 


striction on mobile bed was assumed as 
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(F 5m) (7.6) 


Dimensional analysis of all the pertinent variables also 
Supported the equations 7.4 through 7.6. 

The experimental data indicated the validity of the equations 
7.4 through 7.6. It was found that the backwater rise for the con- 
stricted flow on an alluvial bed is substantially less than that on 
a rigid bed. 

Good correlations between H/h, F and M (Fig. 6.16) and C 
F and m (Fig. 6.17) were noted and design curves have been presented 
to evaluate the maximum backwater depth at any contraction ratio. The 
effect of the variation of the bed material size or the bed scour on 
the maximum backwater rise could not be conclusively studied. It was 
noted that backwater level decreased with increase of bed scour, but 
beyond a certain depth of scour, the backwater level attained a quasi- 
equilibrium condition and no longer decreased with increase in scour 
depth. 

This dissertation synthesizes the procedures for compre- 
hensive hydraulic design of end-dump closure on an alluvial channel. 
Because the design curves agree with prototype data and with the data 
of other researchers, they are considered to be of value for design 


and construction. 
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7.2 Recommendations 
The following additional research studies are considered de- 
Sirable to further the understanding of the complex flow phenomena in 
an end-dump closure and to improve the design procedures: 
(1) Studies on the stability of material on end-dump closure 
dams inclined upstream. It is expected that due to the deflection of 
the high velocity zone away from the advancing face, the closure ma- 
terial would gain increased stability. The total quantity of closure 
material and the bed scour would have to be concurrently assessed. 
(2) Studies to determine the effect of variation of the crest 
width on the stability of closure material. 
(3) Studies to ascertain the effect of rate of dumping more 
decisively. 
(4) Studies to determine the effect of the gradation and shape 
of closure material on the stability. These factors would also influence 
the seepage through the body of the dam, which in turn would change the 
flow through the constriction. 
(5) Detailed studies involving variation of the bed material 
size and gradation over a wider range than heretofore, 
(a) to ascertain the effects on the stability of closure 
material when the bed scour is totally prevented, 
(b) to determine the required size of boulders for paving 


an alluvial bed prior to an end-dump closure, and 
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(c) to enable precise correlations between the pattern 
and extent of scour, and the maximum backwater with 
the bed material characteristics. 
(6) Detailed studies of typical natural river sections to as- 
certain the effects on stability, scour and backwater. 
(7) Since the validity of a design method must finally depend 
upon comparison with prototype performance, and because most river 
closures carried out in the past have not been well documented, there 


is an urgent need for collection and publication of full scale data. 
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APPENDIX A 
CHARACTERISTICS OF NORMAL FLOW 


oi ana Vipers 
ayy a ee . - 7 
4 aha an hy 


ial a pay ryan “oil 


he, 
cf v a f 
i 
PM U n ) 
~ ~ ~ , ii 
: ’ 
if 
an 
— r 


A2 


APPENDIX A 
CHARACTERISTICS OF NORMAL FLOW 


A.1 Introduction 

Simons and Richardson. (1961) have studied and classified the 
bed forms due to flow in alluvial channels into ripples, dunes, anti- 
dunes and plane bed. These bed forms are related to the regime of 
the flow and particle size of the bed material. In general for 
sand-bed channels, after initiation of motion, ripples are formed at 
the low rates of transport. With increasing stream power ripples 
change to dunes and antidunes and eventually back to a plane bed. 
Simons and Richardson observed that "after beginning of motion, 
the plane bed changed to ripples for the sand sizes smaller than 


0.5 mm and to dunes for the 0.93 mm sand". 


A.2 Bed Forms 

Experimental data is precise and collection of data is 
convenient in studies involving clear water flow over a plane bed. 
Theoretical considerations indicated that the effect of variation of 
bed material on the stability of closure material is minor. (Section 
3.1.3). Although a check of the effect of bed material size on 
stability is made in series J through P, the design recommendations 
on stability are basically derived from the data of runs in series A, 


through G. It was decided to conduct these runs with bed material of 
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1.2 mm, for which there was no bed movement under normal flow conditions. 
Design charts on scour and backwater are of course based on the data 

of runs in all series covering three sizes of bed material. With 0.6 mm 
sand, the bed was static and plane ataFroude number of 0.10 whereas at a 
Froude number of 0.29, there was very little movement (plane bed to 
ripples) and at Froude number of 0.50, there was weak bed movement 

(bed form ripples). With 0.25 mm sand, there was weak bed movement 
(ripples at aFroude number of 0.10 and appreciable bed movement (dunes) 
at Froude numbers of 0.29 and 0.50). These observations are similar to 
those of Simons and Richardson (1961) and Cooper (1970). Figure 4.1 
Shows the plane bed of a normal run with 1.2 mm sand. Figure A.1 shows 


the duned bed with 0.25 mm sand for a run in series M. 


A.3 Analysis 


The computed ae: of boundary shear stress ie for the 
normal runs (assuming the entire resistance to be grain resistance 
(Laursen, 1960, 1963) are presented in Table A.1. The computation 
is based on the Manning-Strickler formula which upon algebraic mani- 


pulation yields 
a RF rey (A.1) 
The computation of Gs from Eq. A.1 has obvious limitation 


for a duned or rippled bed, because the total resistance is the result 


of grain resistance and form resistance. The critical tractive stress 
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TABLE A.1 
CHARACTERISTICS OF NORMAL RUNS 


F h V T nc 
Froude Normal Normal Normal Critical 7 //* Bed movement and 
number of depth of velotity Shear Shear Pie form 
approach flow Stress Stress 
flow. in feet in fps (Epc. 1p) 4 (Eqs 3e 29) 


in psf in psf 


Series B through H, d = 1.2 mm 


0.50% 0.18 Wwl5 0.013 0.016 0.785 No movement plane 
bed 

0.29 0.19 0.72 0.005 0.016 Q232 No movement plane 
' 0:26 0.84 0.0065 0.016 0.41 bed 

0.10 0.35 0.38 0.0009 0.016 0.05 No movement plane 
bed 


Series J, K and L, d = 0.60 mm 


0.50 0.18 ier) 0.0097 0.008 eA Slight movement 
ripples 
0.29 0.26 0.84 0.0046 0.008 0.58 Very weak movement 
; plane bed to ripples 
0.10 pie35 O38 0.0007 0.008 0.063 No movement plane 
bed 


Series M, N and P, d = 0.25 mm 


0.50 0.18 1.15 0.0067 0.0033 220) Rapid movement 
dunes 

0.29 0.26 0.84 0.0035 0.0033 1.04 Some movement 
ripples 

0.10 0.35 0.38 0.0005 0.0033 0215 Weak movement 
ripples 


ee 


’The actual water surface slope for the normal flow was 0.0018 against the 
theoretical slope from Manning's formula of 0.00116 for uniform flow. 

(n = 0.013 for d = 1.2 mm). Detailed velocity profiles were recorded for 
the normal flow. From velocities measured at 0.2 depth and 0.8 depth, 

n was calculated assuming a logarithmic velocity distribution from boundary 


layer theory (Chow 1959)as 0.015 
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ie due to Shields (1936) in respect of each material is incorporated 
in Table A.1. It is interesting to see that bed movement was initiated 
at Nae. values as low as 0.15 for 0.25 mm sand. Laursen (1960, 1963) 
also suggests beginning of sediment-transport at Tee values in the 
range of 0.20 to 0.80. 

As indicated on Tables 4.4 and A.1, the runs in series A 
through H and series K were clear water runs and the rest were sedi- 
ment-transporting runs. For the latter runs, sediment was fed manually 
at the upstream end of the flume, but this procedure can not be 
considered to maintain stable approach conditions, as is the case in 
a re-circulating flume. In other words the approach Froude number 
can not be considered to be strictly invariant in the sediment-trans- 
porting runs in this study. This aspect is dealt with in the analysis 
of data in section 6.2.3. It however needs to be mentioned that the 
computed value of the maximum rate of sediment discharge corresponding 
to the Froude number of 0.50 for the bed material of 0.25 mm, according 
to the design charts proposed by Cooper (1970) is only 90 parts per hundred 


thousand, which is indeed a low rate of transport. 
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APPENDIX B 
TURBULENT SEEPAGE THROUGH CLOSURE DAM 
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APPENDIX B 
TURBULENT SEEPAGE THROUGH CLOSURE DAM 


B.1 Introduction 
A closure dam formed of rockfill is a highly permeable 
structure. As closure progresses, the head differential between the 
upstream and downstream of the dam increases, causing considerable 
seepage through the body of the dam. The flow through the large 
voids in a prototype is turbulent. It has been established that 
beyond a particle paynold’s number of about 75, the flow changes 
from laminar to turbulent. Turbulent flow analysis has relied heavily 
on empirical relationships involving the porous material characteristics. 
| In general, the basic turbulent flow relationship has been 
assumed to take the form (Forcheimer, 1930). 
i= ryt ryV + rv" 
(B.1) 


, A 
or 7 = r4V 


in which i is a representative hydraulic gradient and V is the seep- 


age velocity through gross area. 


The evaluation of constants ry> Vo. Ya and ry and the exponent 


A’ has to be done empirically for the particular size and shape of the 


material forming the dam. 
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Taylor (1957) suggested a characteristic linear parameter 
Ra to denote the hydraulic radius of the voids, which is the ratio 


of the volume of voids to the surface area of the particles. 


me volume of voids 
m surface area of particles 


- volume of voids/volume of solids 
surface area of particles/volume of solids 


z porosity era (B.2) 


~ surface area/volume of solids K, 


From studies of Wilkin (1963) and Parkin (1956) the constant Ky 
has been evaluated to be in the range 32 - 48 from tests on crushed 
stones and smooth marbles ranging from 3/8" to 8". Parkin further 


deduced the equation 


in which Ly represents the mean void velocity given by 


Q 
v= Re — (B.4) 
f 


in which Q = seepage discharge, A¢ = gross cross-sectional area 


normal to flow. 


Equation B.3 is similar in form to Chezy's equation 
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i= we (B.5) 


Izbash (1970) analysed all the variables involved in turbu- 
lent seepage from dimensional considerations and derived a similar 
equation 


1/2 
Ve= Cp (y=) [D i] (B.6) 


Ite 
in which Ce, is designated by Izbash as a "generalized Chezy co- 


efficient". 
The equation B.6 is generalized further by Izbash in the form 


1/2 
Vite Kee Lined (B.7) 


in which Ke is the coefficient of turbulent filtration, which is 
primarily a function of D, the nominal size of the material and the 
porosity of fill. It can be seen that equations B.1, B.3 and B.7 
have the same power form and involve constants which require empirical 
evaluation. Relevant studies in this regard have been conducted by 
Hausser and Michael (1959) and Cohen de Lara (1956). 

In Table 5.2, the particularsof runs in which seepage dis- 
charge exceeded 10% of the uncontracted flow, have been presented. 
Assuming that the seepage through the closure dam follows equation 


B.7, an analysis is made herein to determine the coefficient of turbu- 


lent infiltration K- for the closure materials. 
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B.2 Quantity of seepage Through a Closure Dam 


Referring to Fig. B.1, the information available, for the 
nonoverflow portion of the end-dump dam are the head difference, AH, 
the height of the dam,Hp. and the length of the dam at the base,L,. 
Equations B.2, and B.4 or B.7 would not directly yield the seepage 
discharge, as the slope of the phreatic surface is variable in the 
body of the dam and the velocity of flow is not constant either 
longitudinally or vertically. 

Step by step integration of the velocity-slope equation 
will lead to the formulation of total seepage discharge. 

Following Izbash (1970), the saturated zone in the body 
of the dam (Fig. B.1) can be considered to be made of two parts, a trape- 
zium ABDE belowa horizontal line at the downstream water level and 
a triangle BCD above it. 

From Dupuit's theory, which basically assumes that for small 
inclinations of the line of seepage the streamlines can be taken as 
horizontal, (Harr, 1962), the flow velocity is given by (Eq. B.7) 

= Ke dy (B.8) 

By integration of Eq. B.8 with respect to depth the total 


seepage through the trapezoidal part is given by 


H, AH , 
ag = Bhp YE [My - Ay = hy] (B.9) 
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and seepage through the triangular part is given by 
s (AH 
te Kev 31Hh (B.10) 


te is generally much smaller than Ve - Total seepage per unit width 
] 


of the dam is given by 
eee (B.11) 


B.3 Computation of Ke and Discussion 


The hydraulic particulars(in equations B.9 and B.10) for the 
runs in Table 5.2 are presented in Table B.1. (Fig. B.la is to be 
referred as the key diagram.) From the known values of qr for each 
run, Ke is computed and presented in Table B.1. 

A plot of Ge versus 2 for the three sizes of closure ma- 
terial is shown in Fig. B.1b. Inspiteof the scatter, the plot indicates 
a definite correlation of q¢ and z. The slope of the line gives Ke. 

As shown by Parkin (Eq. B.3), Ke does not really have a unique re- 
lationship with the particle size only and depends upon porosity and 
shape factors which are possible reasons for scatter of the data in 
Rigo. 1b: 

The mean values of Ke computed from the data in Table B.1 works 
out to 0.13, 0.36 and 0.56 ft/sec for the three closure materials of 


size 0.021, 0.059 and 0.082 ft. respectively. A logarithmic plot of 
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Ke versus D shows a trend (Fig. B.2). Izbash's values of Ke for stone 
chippings are plotted on Fig. B.2. Figure B.2 is not proposed to be 
extrapolated or used as a design curve. 

It appears that a power form of relationship between seepage 
velocity and the hydraulic gradient (such as equation B.7) works 
reasonably in the turbulent flow regime over the size range tested 


in this study. 
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APPENDIX C 
ANALYSIS OF DATA BY SHEAR-STRESS STABILITY CRITERIA 


C.1 Introduction 

As observed in this study (section 5.5.1), the stable size 
of the closure material did not agree with that given by correlation 
(eq. 2.4) with either the mean velocity through the constriction or 
the maximum velocity against the closure material. It has been esta- 
blished in section 5.5.2 that the experimentally observed stable size 
agrees satisfactorily (Fig. 5.17) with the theoretical analysis, 
(section 3.1.3), in which the basic assumptions were 

(i) stable size can be determined by correlation with the 
maximum drag = on the face. 


max 


(ii) ies [te increases exponentially with the increase of 
max 


contraction ratio. 

(Shear stress and drag are synonymous, denoting the hydrodynamic drag 
impaistinction to litt). 

In the course of this investigation detailed velocity 
measurements normal to the end-dump face have been recorded from 
which drag distribution on the face are obtained by assuming a 
boundary ees flow. Analysis and examination of the data at thres- 


hold condition are attempted herein from shear-stress stability 


criteria. 


G2 


. — 
' , 245 i od 
ALASTIAD VI LIIBATS. 4 


psy 


~ 


este wldate one ft ed,2 not rase) tbe tae att saad 
notisiovims yd favre sary ttiw aertye J00 ‘pte: fetaseen, 
40 nots .tano! etd feuords woofs fag _old sossia data 
~st29 nad asd TT etaedae oweofs ode rentsps wWiretey ax 
exhe etdata bawtardo \ Karmen veqxe ey todd $22 wn 8 
Lcteyfone! Tpohipyaamt alg Ashe (4-8 at) ah 

9 thw snot tqmvess. ahead seta sate, 

oft dttw notislories hawtnvatob ad asa ssto-otdate 
cao? efit ao at pate a 7 

io sessvoer ont dttw et ferdwswaqne: anneal - - 


pat otnangborinye: alt aiisondlb <2uemgoome e18 arb ban sande 
| CAT 08 notzonttatt itt 

Tootey ber hadsb moktepEsement ata Yo eenyoo aly At 4 

oo? bebro981 asad\sved sos gub-tnd eit? of Teorion ednamauaiem 

< ontaweas vt bentasdo ove sont ot na malgudtayats web doi” 

220i? 26 ateh st Yo eotsehhmexs bis: dheytaiiA wot Tayat ysbmued 

ustlidete eesvde~teariz mort ate wil botgmaiss ev6 mote tbnos Blot | 


C3 


C.2 Experimental Measurement of Drag on Advancing Boundary 


The shear stress in the flow contraction zone can be deter- 
mined using the Karman-Prandtl velocity distribution for fully 


developed boundary layer flow over a rough bed, given by 


{= 5.75 log jee 8.5 (oni) 
~ S 


in which u represents the velocity at a height y above the surface 
and u, is the shear velocity and Kes the roughness. 

The constant 5.75 is generally applicable and constant 8.5 
can be assumed for turbulent flow with particle Reynolds Number above 
90, which was the case with the coarse closure material. The appli- 
cation of boundary layer theory is justified because the flow enters 
the contraction under a favourable pressure gradient and does not 
leave the boundary in the zone under investigation. The velocity dis- 
tribution normal to the curved boundary can therefore be assumed to 
follow Equation C.1 . Eliminating k, and constant 8.5 by algebraic 
manipulation, drag on the boundary can be evaluated from the equation 


(ae 


5.75 109 Yo/Vy (C.2) 


where U, and Uy represent flow velocity at depth Yo and yy normal to 


the boundary respectively. Figure C.1 shows the definition sketch and 


a typical velocity-depth semilog plot. 
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By resorting to the above procedure the shear stress distri- 
bution in the contraction zone as well as the maximum drag on the end- 
dump face Poa) were obtained. A typical plot of non-dimensionalized 
shear distribution (contours of eat) is shown in Fig. C.2. In 
section A, ues the shear stress on the bed for uncontracted or normal 
flow has been evaluated for each run in this study. 

Figure C.2 shows that shear stresses many times the normal 
occur on the advancing boundary. The locus of maximum shear has a 
shape similar to the live stream boundary as does the maximum velocity 
filament. The shear on the backwater reach and in the bed at the 
dump line is well below that on the slope of the dam. All these features 
parallel Ippen and Drinker's (1962) observations for flow along a 
trapezoidal channel b nd. 

The shear stress maxima were computed for each run. A semi- 
log plot of non-dimensionalised shear maxima on the advancing boundary 
versus the contraction ratio is shown in Fig. C.3 for all the runs in 
series B, C and D. The considerable scatter of the data on Fig. 3053 
is to be expected, because the shear measurement based on the assump- 
tion of boundary layer flow over rough, sloping and curved boundaries, 
has obvious limitations. 

Figure C.3, however, clearly indicates that 

(i) In the range tested, drag on the boundary increases ex- 


ponentially with the contraction ratio. This attests. to 


the validity of equation 3.6. 
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NOTE: IN THE SHADED REGION FLOW 
SEPARATED FROM THE BOUNDARY 


FIG. C.2 NON-DIMENSIONAL BOUNDARY SHEAR-STRESS 
(cof) DISTRIBUTION IN CONSTRICTION AT m = 0.56 
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FIG. C.3 VARIATION OF NON-DIMENSIONALISED SHEAR MAXIMA 
IT: ON ADVANCING BOUNDARY WITH CONTRACTION RATIO m 
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(ii) The effect of the variation of Froude number of approach 


flow on t /t. is not significant. 
Cae aie 


As T) is known from the approach flow data, the mean curve 
in Fig. C.3 is used to determine 18 approximately at any stage of 


max 
closure. 


C.3 Effect of Hydrodynamic Lift 
A.S.C.E. Task Committee on Sedimentation (1966) state "most 


works on initiation of motion consider only shear stress and completely 
ignore lift on particles, despite the fact that lift must be present”. 
Based on the studies of Einstein and El-Samni (1949) the Task Committee 
concluded "that the lift is of considerable importance in entraining 
sediment". Studies of Chepil (1958) substantiates this. The essential 
feature of the studies is that lift on the top layer of particles, like 
drag, is also proportional to the square of the velocity close to the 
particle. (Einstein found U5 which is the velocity measured at 0.35 

igs above theoretical bed as representative). 

Recent studies of Benedict and Christensen (1971), and 
Christensen (1971) indicate that flows in the hydrodynamically rough 
range will have positive values of lift. Christensen (1971), studying 
the state of incipient motion on cohesionless channel banks argues 
that "prediction of instability using the classic USBR equation 
based on critical drag criteria (Carter, 1953, and Lane et al 1953) and 


neglecting lift, must be expected to be somewhat on the unsafe side". 
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C9 
Christensen developed a new law for the velocity profile 
close to the bed (where logarithmic profile breaks down by yielding 
negative time-mean velocities) and deduced therefrom a relation for 
the lift-shear stress ratio. Christensen ultimately proposes a 


complex form of modification to the bank correction factor (equation 3.8) 


as 
10.4D 2 
(s-cot“s) (coto+0.556[1(=— #1)] ) 
“eb(lift) - 
1g 
C 10.4D 2 10.4D 4 
0.556s[ In(——> +1) ] +v(s¢-cot@)+0.309[1n a= +1) 
S S 
| (c.3) 


cotoY 1ts¢ 


in which s = slope of the bank. 

This correction factor is always less than the value ob- 
tained from equation 3.8. It will, however, approach that value 
when k./D3, approaches infinity. Equation C.3 has been given in 
graphical form for easy design use. Christensen's recommendation 
being relevant to the situation on the sloping end-dump face, is applied 


herein to determine the lift on the top layer of particles. 


C.4 Stability Analysis Considering Drag and Lift 


To test the hypothesis that the closure material gets en- 


trained corresponding to a critical combination of the drag and lift 
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forces, the following computational steps are followed. 

(1) From the design chart in Fig. 5.15, contraction ratio cor- 
responding to a closure efficiency of 0.90 (acceptable stable condition 
as recommended in this study) is read out for each size of closure 
material at each of the three Froude numbers of 0.50, 0.29 and 0.10. 

(2) Sah corresponding to the contraction ratio in Step 1 is 
obtained from Fig. C.3, using the mown value of Toe (Table A.1). 


(3) te the critical shear stress for the closure material on 


a horizontal bed is computed from 
oes 4D in fps units (E. a) 


(4) As recommended by Christensen, T the actual shear 


max 


stress occurring on the face should be equal to Tob (lift): oe 


is calculated for each case. 

(5) As $, the angle of repose of each closure material, and s, 
the slope ~ 2(horizontal): 1(vertical) is known, the design curves 
of Christensen based on equation C.3 are used to obtain the critical 
size of closure material corresponding to oe (6 of the critical 
size obtained, should strictly be used and not 9 of the actual closure 
material in use of the design curves. The error introduced is how- 
ever insignificant). 


(6) The critical size Dy, obtained is compared with the actual 


size D, corresponding to design stability condition. 
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The observed and computed data following the above steps is 
presented in Table C.1. In Table C.1, the critical size D,, corres- 
ponding to Poel (neglecting lift) based on the USBR equation is also 
presented in Column 8. 

The agreement between D, and D is indeed satisfactory, 
except at m values exceeding 0.80 (Series C). This is due to the 
possible deviation from the exponential relation between To ay 10 
and m, after the two wings of the dam join at the centre which results 
in abnormal increase of both velocity through the constriction and 
Shear stress on the face. 

The inadequacy of neglecting lift in sediment entrainment 


‘is obvious from the consistent disagreement between D3(Col. 8) and 


Colt |): 


C.5 Analysis by Shields' Dimensionless Shear Stress 


Shields (1936) experimental investigations of flow on a 
level bed, led to the generally accepted entrainment plot of the non-di- 
mensionless shear stress Ty = To/(9,-P,,)D : versus boundary Reynolds 
number Rx = uxd/v d For turbulent flow (beyond a Reynolds number 
of 1000) t, appears to be independent of R, and takes a value of 0.04 
to 0.06. 

Although application of Shield's diagram to the incipient 
condition on the curved, sloping face of the end-dump is not justified, 
an analysis is made by plotting Toy! PD Pw? as abcissa against 


efficiency of closure, n as ordinate in Fig. C.4. The data plotted 
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C14 
pertain to the runs in Series B, C and D, for which the efficiency of 
closure was above 0.70. (The measurement of t at closure efficiencies 
below 0.70 is considered erratic, because of A VeTTOnE of material on the 
boundary). 

It is interesting to note that even for highly curvilinear 


approach flow, Shields parameter takes a value of 0.02 at an efficiency 


of 0.90, considered as stable design condition. 


C.6 Concluding Remarks 


Maximum shear stressesoccurring on the boundary have been com- 
puted from measurement of velocity normal to the end-dump face. Good 
agreement is observed between the critical size and the value computed 
by Christensen's curves correlating the size with drag and lift. An 
analysis in the light of Shields diagram indicated t, to take a value 
of 0.02 for the curvilinear flow on end-dump face. Correlation of 
stable size with drag (neglecting lift) was found to be unsatisfactory. 
Correlation of stable size with the value given by critical velocity 
criteria was observed to be unrealistic (section 5.5.1). 

Although the above observations are not conclusive regarding 
the practicability of selecting material size on the basis of shear 


stress, they convincingly indicate abandoment of the velocity criterion 


in favour of the proposed approach. 
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